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ABSTRACT 


This  thesis  deals  primarily  with  the  development  of 
an  approximate  method  of  analysis  of  a  tunnel  diode  oscillator. 
Because  of  the  non-linear  voltage-current  characteristic  of 
the  tunnel  diode,  relatively  simple  systems  may  exhibit  cer¬ 
tain  phenomena  which  cannot  be  explained  in  terms  of  linear 
circuit  theory. 

If  two  mutually  coupled  parallel  resonant  circuits 
comprise  the  linear  portion  of  a  tunnel  diode  oscillator, 
under  certain  conditions,  two  distinct  modes  of  oscillation 
will  be  exhibited.  Furthermore,  as  the  degree  of  coupling  is 
varied,  two  critical  points  can  be  found  for  which  the  oscil¬ 
lator  will  switch  from  one  mode  to  the  other.  A  form  of 
hysteresis  is  also  exhibited  in  that  the  point  at  which  the 
oscillator  will  switch  to  its  second  mode  occurs  for  a  larger 
degree  of  coupling  than  does  the  point  at  which  it  will  switch 
back  to  its  original  mode. 

The  objective  of  this  thesis  is  thus  to  present  a 
method  of  analysis  sufficiently  accurate  as  to  account  for 
the  above  phenomenon.  This  is  done  by  first  analysing  an 
oscillator  using  a  single  tank  circuit,  and  then  considering 
the  effect  of  the  second  tank  circuit. 
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I .  INTRODUCTION 
1 . 1  The  Tunnel  Diode 

In  the  analysis  of  any  non-linear  system,  the  fac¬ 
tors  of  primary  importance  are  the  characteristics  of  the  non¬ 
linear  device.  Thus,  only  a  brief  qualitative  discussion  of 
the  operation  of  the  tunnel  diode  itself  will  be  presented.  A 
more  sophisticated  treatment  can  be  found  in  a  variety  of 
literature . 

The  tunnel  diode  is  a  semiconductor  device  contain¬ 
ing  a  single  p-n  junction.  It  differs  primarily  from  the  con¬ 
ventional  semiconductor  diode  in  that  the  levels  of  doping  in 
both  the  p-type  and  n-type  materials  are  many  times  higher. 

This  high  concentration  of  impurities  results  in  an  extremely 
thin  depletion  region  at  the  junction.  For  very  low  values 
of  bias,  this  depletion  region  represents  an  energy  barrier 
which  electrons  are  unable  to  surmount.  However,  because  of 
the  tunnel  effect,  whereby  a  particle  may  disappear  from  one 
side  of  such  a  barrier,  and  appear  on  the  other  side  without 
loss  of  energy,  an  appreciable  current  will  flow  across  the 
junction.  As  bias  is  subsequently  increased,  the  tunnelling 
current  increases  up  to  a  point,  and  then  decreases  to  approxi¬ 
mately  zero.  This  gives  rise  to  the  "negative  resistance' 
region  of  the  tunnel  diode.  Upon  further  increase  in  bias, the 
electrons  are  given  sufficient  energy  to  surmount  the  barrier. 
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and  conduction  takes  place  by  means  of  the  same  mechanism  as 
in  an  ordinary  diode. 

1 . 2  Tunnel  Diode  Characteristic 

The  d.c.  voltage-current  characteristic  for  a  typi¬ 
cal  tunnel  diode  is  illustrated  in  Fig.  1.1.  The  meanings  of 
a  number  of  important  parameters  are  also  indicated. 


Figure  1.1 

Tunnel  Diode  Characteristic 

The  functional  notation  I  =  f (V)  shall  be  used  for 
a  mathematical  description  of  the  d.c.  characteristic.  With¬ 
in  the  region  between  V  and  V  ,  where  f ' (V)  is  negative,  the 
^  p  v 

tunnel  diode  can  be  thought  of  as  a  small-signal  negative  re¬ 
sistance.  Hence  it  can  be  used  as  a  source  of  a.c.  power,  and 
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readily  lends  itself  to  amplifier  and  oscillator  applications. 

The  quantity  r^  shall  be  defined,  in  the  interval 
<  V  <  Vv,  by  the  equation 

1 

-  —  =  f 1 (V)  .  (1.1) 


Also,  the  quantity  -G^  shall  be  defined  as  the  minimum  value 
of  f ' (V) ,  or 


-Gd  =  f,(VI)  .  (1.2) 

The  symbol  to  be  used  throughout  this  thesis  for  a 
tunnel  diode  is  shown  in  Fig.  1.2.  The  normal  biasing  polar¬ 
ity  is  also  indicated. 


o 

Figure  1.2 

Symbol  for  a  Tunnel  Diode 

1 . 3  Parasitic  Elements 

Any  tunnel  diode,  when  used  in  conjunction  with  ex¬ 
ternal  circuitry,  will  have  associated  with  it  parasitic  cir¬ 
cuit  elements,  R  ,  L  ,  and  C  ,  where: 

'  s  s  s 

Rg  represents  a  series  resistance  due  to  the  bulk 
material  of  the  diode  and  associated  wiring; 
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Lg  represents  a  series  inductance  mainly  due  to 
the  wiring?  and 

C  represents  a  shunt  capacitance  depending  on  the 
s 

junction,  the  package,  and  the  wiring. 

If  an  "ideal”  tunnel  diode  is  now  defined  as  a  de¬ 


vice  which  produces  only  the  d,c.  characteristic,  I  =  f (V) , 
the  parasitic  elements  can  be  lumped  and  connected  externally 
to  this  ideal  diode.  The  resulting  network  can  be  used  as 
the  equivalent  circuit  for  a  "real"  tunnel  diode.  This  is 
shown  in  Fig.  1.3. 


C 


S 


Figure  1.3 

Equivalent  Circuit  for  a  Real  Tunnel  Diode 

Note  that  since  the  series  resistance  associated 


with  the  tunnel  diode  itself  will  be  included  in  determining 
the  d.c.  characteristic/  Rg  /  as  indicated  in  Fig.  1.3,  repre¬ 
sents  only  the  series  resistance  external  to  the  actual  diode. 


In  the  analysis  of  the  tunnel  diode  oscillator,  the 


parasitic  elements  have  first  been  neglected,  and  their  pos¬ 
sible  consequence  are  later  considered  in  Sec.  5.10. 
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1.4  Basic  Oscillator  Circuit 


Figure  1.4  shows  the  basic  form  of  oscillator  to  be 

analysed  in  this  thesis.  The  biasing  network  (V__ ,  R_  ,  R_  , 

BB  B  B  0 
1  £ 

and  Cg)  is  used  to  bias  the  tunnel  diode  in  the  'negative  re¬ 
sistance*'  region.  It  will  be  assumed  that  Cg  is  sufficiently 
large  that  its  impedance  can  be  neglected  at  the  frequency  of 
oscillation . 


Figure  1.4 

Basic  Oscillator  Circuit 


If  the  secondary  tank  circuit  (Fig.  1.4(a))  is  first 
neglected,  the  operation  of  the  oscillator  can  be  qualitatively 
understood  by  considering  the  tunnel  diode  as  a  source  of  a.c. 
energy,  which  compensates  for  the  energy  losses  in  the  primary 
tank  circuit. 

When  the  secondary  tank  circuit  is  mutually  coupled 
to  the  primary,  the  coefficient  of  coupling  can  be  defined  as 


K 


2 


M 


2 


L1L2 


/ 


(1.3) 
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where  M  is  the  mutual  inductance  between  and  .  Thus, 
varying  the  distance  between  and  L 2  will  result  in  an 
opposite  variation  in  . 

The  nature  of  the  switching  phenomenon,  as  discussed 
in  this  thesis,  is  qualitatively  illustrated  in  Fig.  1.5.  In 
Fig.  1.5(a),  f  refers  to  the  fundamental  frequency  of  oscilla¬ 
tion,  and  in 


Fig.  1.5(b),  A  refers  to  the  fundamental  amplitude  of  time- 

varying  voltage  across  the  tunnel  diode. 

As  can  be  seen  from  Fig.  1.5,  two  distinct  modes 

of  oscillation  are  exhibited,  and  there  are  two  critical 

values  of  K  for  which  the  oscillator  will  switch  modes.  If 

K  is  between  these  two  values,  both  modes  are  available,  and 

2 

the  mode  assumed  by  the  oscillator  depends  on  whether  K  has 
recently  been  increased  from  a  small  value  or  decreased  from  a 
large  value. 
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A  further  experimental  observation  is  that  the  pat¬ 
terns  shown  in  Fig.  1.5  seem  to  be  exhibited  only  when  the 
secondary  tank  circuit  has  a  considerably  higher  Q  and  a 
slightly  lower  resonant  frequency  than  the  primary. 

1 . 5  Methods  of  Analysis 

Two  methods  of  analysis  of  the  oscillator  circuit 
are  presented.  The  first  method  considers  the  tunnel  diode 
oscillator  as  a  feedback  oscillator.  This  method  essentially 
linearizes  the  d.c.  characteristic  at  the  point  of  operation, 
and  is  thus  valid  only  for  very  small  amplitudes.  It  can  be 
used,  however,  to  derive  certain  necessary  conditions  for 
oscillation . 

The  second  method,  which  is  much  more  fruitful  than 
the  first,  is  essentially  the  describing-function  method. 

The  basic  assumption  of  this  method  is  that  time-varying  volt¬ 
age  across  the  tunnel  diode  is  purely  sinusoidal.  The  essen¬ 
tial  differences  between  this  method  and  the  conventional 
describing-function  method,  are  that  (1)  the  describing-func¬ 
tion  equation  is  a  function  of  both  amplitude  and  operating 
point,  and  (2)  the  shift  in  operating  point  due  to  a  change 
in  amplitude  is  considered. 
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II.  FEEDBACK  OSCILLATOR  MODEL 
2 . 1  Equivalent  Biasing  Network 

If  the  d.c.  bias  voltage,  VDtl,  in  Fig.  1.4(a),  is 

JdJd 

applied  suddenly  at  some  instant,  it  can  be  considered  as  a 
voltage  step,  which  has  the  Laplace  Transform  V_._./S .  The 

DO 

biasing  network  can  thus  be  represented  as  shown  in  Fig. 2. 1(a) 


R 


B. 


VS) 


R 


B 


WWV 


'B 


(b) 


Figure  2.1 

Equivalent  Biasing  Network 

Application  of  Thevenin's  Theorem  now  results  in 
the  equivalent  circuit  shown  in  Fig.  2.1(b),  if 


R 


B 


Vb2 

Rb1  +  **2 


(2.1) 


Simple  analysis  shows  that 


VBB  RB. 


vs)  = 


S(R„  +  Rn  +  SC  R  Rr  ) 
B1  B2  B  B1  2 


(2.2) 
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If  the  analysis  is  now  confined  to  steady-state 
conditions,  Vp(S)  in  Fig.  2.1(b),  can  be  replaced  by  its 
steady-state  value.  Using  the  Final  Value  Theorem  of  Laplace 
Transforms , 


V 


B 


lim  VB(t) 

t +°° 

vbbV 


lim  SVD (S) 

S-0  B 


rr  +  rr 
B1  B2 


(2.3) 


Thus,  the  circuit  of  Fig.  1.4(a)  can  be  replaced  by  it$ 
"steady-state1'  equivalent,  shown  in  Fig.  2.2(a). 


R. 


+  VD  - 


Figure  2.2 

Steady-State  Equivalent  Circuits 


The  circuit  of  Fig.  2.2(a)  can  in  turn  be  represented 
by  Fig.  2.2(b),  if  ZT  represents  the  tank  circuit  impedance, 
and  ZD  the  biasing  network  impedance. 

±3 


Si. 


vd  3d  neo  <e)*.I  .pi1?  5o  diuouio  arti  ,a«ri' 

>  '  ‘  1  •' 
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2 . 2  Small  Signal  Equivalent  Circuit 


For  very  small  signals  about  some  operating  point, 
the  tunnel  diode  can  be  approximated  by  a  linear  negative  re¬ 
sistance,  -r^.  Furthermore,  it  has  already  been  assumed  that 

at  all  frequencies  of  interest,  Z  =  0.  Thus,  assuming  that 

a 

all  d.c.  currents  and  voltages  remain  constant,  the  circuit  of 
Fig.  2.2(b)  can  be  replaced  by  its  small  signal  equivalent, 
shown  in  Fig.  2.3(a). 


+ 


Figure  2.3 

Small-Signal  Equivalent  Circuits 


Note  that  for  generality,  an  external  signal,  v^ ,  has  been 
added . 

Two  equations  can  be  written  for  the  circuit  of 
Fig .  2.3(a): 


v  ^  —  i^Z(p  -t-  v^  r  and 


(2.4) 


(2.5) 


y 
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Referring  now  to  Fig.  2.3(b),  if  6  =  -1,  the  same  two  equa¬ 
tions  apply,  so  that  this  circuit  will  be  identical  to  that 
of  Fig .  2.3(a). 

The  circuit  of  Fig.  2.3(b)  is  easily  recognized  as 
a  feedback  amplifier,  with  a  feedback  factor  3,  and  an  open- 


v 


loop  gain  A 


T 


v 


v. 


v. 


For  B  =  0,  i,  =  — 1 

1  -r 


and  vT  =  i1ZT. 


3=0 


Hence , 


A  =  - 
v 


(2.6) 


Now  setting  3  =  -1  is  essentially  closing  the  feedback  loop, 
and  the  closed-loop  gain  of  Fig.  2.3(b)  will  be  given  by 


A, 


v 


1  -  3A 

v 


ZT^rd 

1  -  Vrd 


(2.7) 


This  must  also  be  the  gain  of  the  circuit  in  Fig.  2.3(a), 

since  the  two  circuits  are  identical  for  3  =  -1. 

Since  for  oscillation  of  a  feedback  amplifier, 

3A  =  1,  it  follows  that  the  circuit  shown  in  Fig.  2.3(a)  will 
v 

oscillate  if  ZT/r^  =1,  or 


=  r 


(2.8) 


If  oo  is  now  defined  such  that  Im 
o 


ZT (3“0) 


=  0, 


and  Rq  is  defined  as  Rq  =  ZT(jcoQ),Eq.  2.8  implies  that  oscil¬ 
lation  of  frequency  ooQ  will  occur  if 


R 


=  r 


(2.9) 


. 


n  (,,.10  nM;  c  -baaolo  ^ 


■ 
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Note  that  if  Rq  is  greater  than  r^,  the  open-loop 
gain  will  be  greater  than  one,  and  upon  closing  the  loop,  oscil¬ 
lations  will  build  up  until  some  form  of  limiting  occurs.  This 
is  what  normally  occurs  in  a  tunnel  diode  oscillator. 
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III.  PROPERTIES  OF  THE  PARALLEL  RESONANT  CIRCUIT 
3.1  Resonance 


Since  the  oscillator  can  be  expected  to  oscillate 

with  frequency  ooQ,  such  that  ZT(ja)Q)  =  ,  this  frequency 

shall  be  referred  to  as  the  resonant  frequency  of  the  tank 

circuit.  Similarly,  Rq  shall  be  referred  to  as  the  resonant 

resistance.  Note  that  for  this  circuit,  w  does  not  correspond 

o 

to  an  extremum  of  |  ( j  cd )  |. 

The  impedance  of  the  primary  tank  circuit  can  be 
expressed  in  terms  of  S  as 

R.  +  SL, 

Z  (S)  =  - - - - -  ;  (3.1) 

S  L  C.  +  sciRi  +  1 

or  in  terms  of  jo)  as 


R 


Zm(ja))  = 


2  + 


(1  —  0)  L-^C-^) 


(3.2) 


r  i 


Setting  Im 


Zm  ( i  a)  ) 
T  J  o 


=  0,  it  can  easily  be  shown  that 


R. 


0) 


Lici 


(3.3) 


Substituting  (3.3)  into  (3.2); 


R 


Rici 


(3.4) 


. 


.|{ut)TS|  3o  munisT^xs  rrs  oi 


(2) 
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3 . 2  Frequency-Normalized  Impedance 


Normalized  frequency,  u,  shall  be  defined  such  that 


2  2 
u  =  oo  L^C^ 


(3.5) 


The  tank  circuit  impedance  can  now  be  expressed  in  terms  of 
ju  as 


Zrn(ju) 


Rl(1  +  j  Rl  *VC1> 


(1  -  u  )  +  ju 


(3.6) 


If  Q-^  is  defined  as 


1 

R 


1  V  CX 


(3.7) 


Eq.  (3.6)  can  be  rewritten  as 


ZT(  ju) 


R^  ( 1  +  jQ-^u) 

“  27  ^  .  u 

(1  -  u  >  +  3  07 


(3.8) 


With  the  above  definition  of  in  mind,  (3.3)  can  be  re¬ 
written  as 


0) 


/Lici 


rl  - 


(3.9) 


Q- 


from  which  the  normalized  resonant  frequency  will  be 


»  .  AT2 

°  /  Q, 


(3.10) 


as  ndddiiwr*!  sd  nso  (d.^  .|  i 


l 
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Thus,  if 


Qi  >>1 


oj 


o 


1 

/l7c7 
±  1 


and  u 

o 


1. 


3 . 3  Frequency  Response  Above  Resonance 


Equation  3.8  can  be  used  to  determine  approximately 
the  impedance  to  any  harmonic  of  the  resonant  frequency.  The 
square  of  the  magnitude  of  (3.8)  will  be  given  by: 


ZT ( ju) 


2  2  2 
R1  (1  +  Q1  u  ) 

,,  2^  2  ^  u2 
(1  -  u  )  +  JT 


(3.11) 


If  it  is  assumed  that  Qn  >>1,  u  -  1.  Thus,  for 

1  o 

th  > 

the  n1"  harmonic  above  resonance,  u  -  n.  Since  n  >_  2 ,  u  ~  2, 

2 

2  2  2  2  m  ^ 

and  thus  Q,  u  >>1.  Also,  (1  -  uz)  >>— j  .  Hence, 


2^  2  2 

i 2  i  i 2  R1  ^1  n 

ZT(jnuo)  |  -  |ZT(jn)  |  -  2  ,.2 

(n  -1) 


n 


2  2  ' 
C1  (n^  -  l)z 


or 


ZT ( jn) 


n 


Cx  (n  -  1) 


(3.12) 


th 

Thus,  the  magnitude  of  impedance  to  the  n  harmonic  above 
resonance  is  nearly  independent  of  R-^  (if  Q-^  is  high)  ,  and 
depends  almost  directly  on  /L^/C^  .  The  significance  of  this 
will  be  seen  later. 

Also  the  ratio  of  the  magnitude  of  the  impedance  to 
the  n*"  harmonic  to  the  impedance  at  resonance  can  be  found; 
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or 


But 


|  Z ( jnuQ)  1 

|  Z ( jn)  | 

/h 

/ 

|2(juo>  I 

R 

o 

/  C,  (n2  -  1)  R 

1  o 

1  A 

R  /  C, 

o  1 

i 

Qi  ' 

so  that 

|z(jnuQ) 

n 

|Z<juo) 1 

(n2-l)Q1 

(3.13) 


3 . 4  Bode  Plot 

An  approximate  Bode  Plot  of  |ZT(ju) |  will  now  be 
derived.  This  is  shown  in  Fig.  3.1.  Note  that  absolute  mag¬ 
nitudes  of  impedance  (in  db)  are  used,  so  that  | ZT ( j u)  |  =  1 
corresponds  to  0  -  db.  Equation  3.11  can  be  used  to  determine 
the  as^ymptotes,  critical  frequencies,  and  magnitudes.  The 
following  steps  are  required  to  derive  the  Bode  Plot: 

(1)  If  u2<<1,  |  ZT  ( ju)  |  2  -  R-l2(1  +  Q12u2)  .  This 

has  a  corner  frequency  at  u  =  1/Q^ . 

For  u  <  1/Q^,  the  as^ymptote  is  a  horizontal 
line  of  ordinate  value  20  log^  . 

For  u  >  1/Q^,  the  as^ymptote  will  be  deter¬ 
mined  by  the  equation  |ZT(ju)  |  -  R-^Q^u  =  u  /L^/C^.  This  as- 

y§ymptote  will  have  a  slope  of  +20db/decade ,  and  will  inter¬ 
sect  the  0  -  db  line  when  u  /L^/C^  =  1,  or  u  =  /C^/L^  . 


. 
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R  Q 

(2)  If  u^>>1,  |z  (ju)  |  -  -----  =  —  /l,  /C-,  .  The  cor- 

T  u  u  V  1 

responding  as/lymptote  will  have  a  slope  of  -20db/decade ,  and 
will  intersect  the  0  -  db  line  at  u  =  /Cj  . 

(3)  The  as^ymptotes  having  positive  and  negative 
slopes  will  intersect  where  u/L^/C^  =  ~  /L^/C^  ,  or  u  =  1 . 

At  this  point  the  ordinate  value  (on  either  as^ymptote)  will 

be  20  ^Li/Ci  =  10  log^g  Ll//Cl  * 

(4)  At  u  =  1 ,  |  Z  ( j u )  |  -  R,Q,2  =  R  .  Thus  the 

1  XX  o 

ordinate  value  of  the  actual  curve,  at  u  =  1,  will  be  approxi¬ 
mately  20  log^Q  Rq  =  20  log^Q  L^/R^C-^.  The  approximate  Bode 
Plot  can  now  be  sketched. 


Figure  3.1 

Approximate  Bode  Plot 


. 


X— «  |  U»t)„S|  ,I«  u  51  *S) 


i-  n  »V,  :  'V  *i\  =* 
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Figure  3.1  can  be  used  to  study  the  effect  of  vary¬ 
ing  the  circuit  components.  It  is  first  noticed  that  the  only 
two  separate  factors  appearing  in  the  derivation  of  the  Bode 
Plot  are  and  (L^/C-^)  .  (i.e.  and 

Ro  -  SJ 

If  it  is  desired  to  increase  this  can  be  done 

by  either  decreasing  or  increasing  (L^/C^) .  If  R^  is  first 
decreased,  it  is  apparent  from  Fig.  3.1  that  the  horizontal 
as^ymptote  will  be  shifted  downward,  but  the  two  sloping 
as/ymptotes  will  remain  fixed.  Furthermore,  the  peak- point 
will  be  shifted  upward.  Thus,  near  the  resonant  frequency, 
impedance  will  be  increased;  for  low  frequencies,  impedance 
will  be  decreased;  and  for  high  frequencies,  impedance  will 
remain  nearly  fixed. 

If  R^  is  now  held  constant,  and  is  increased  by 
increasing  L^/C^,  the  horizontal  as^ymptote  will  remain  fixed, 
and  the  sloping  as^ymptotes  will  be  spread  apart.  Thus,  the 
ordinate  of  intersection  of  these  as^ymptotes  will  be  raised, 
and  the  peak-point  will  correspondingly  be  raised  by  approxi¬ 
mately  twice  as  much  (in  db) .  However,  all  high-frequency 
impedances  will  also  be  raised  by  approximately  half  as  much 
(in  db)  as  is  the  peak-point.  It  is  shown,  in  Sec.  5.11,  that 
increasing  in  this  manner  may  result  in  an  increase  in 
harmonic  distortion,  rather  than  a  decrease. 


. 


:i  b  *  d  rep  ai.  ^  ,X'J  easenoai  od  si  di 


/<  >i  9d  v.  ?*  b  ■•.  I  d;  t<  .**  ^<3  ®  " 

■*. 

sioni  ne  ni  dl-x/adi  yem  lertnj?  axirid  nx  pnx  x  ^ 

’ 

.sBB9?03b  s  nsrli  isrfdB^  ,noxd^odexb  oiflounted 
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IV.  DESCRIBING-FUNCTION  METHOD 
4 . 1  Quantitative  Description  of  the  D.C.  Characteristic 


The  analysis  to  follow  is  numerically  based  on  a 
1N2940  Germanian  tunnel  diode.  Table  4.1  lists  the  typical 
parameters  for  this  diode.  These  values  have  been  taken  from 
the  first  edition  of  the  "General  Electric  Tunnel  Diode 
Manual".  (Ref. 3) 


Table  4.1 


Typical  Parameters 


I 

1.00 

ma 

p 

I 

0.13 

ma 

V 

V 

60.0 

mv 

P 

V 

350.0 

mv 

V 

vfP 

500.0 

mv 

Gd 

6.60 

mmho 

Before 

beginning  any 

sort  of 

mathematical  description  of  the  curve  I  =  f (V)  must  be  found. 

t  h 

Such  a  curve  can  be  approximated  by  a  K  order  polynomial  .in 
V,  where  the  value  of  K  as  chosen  depends  on  the  degree  of 
accuracy  required, 
i .  e . 


I  =  f (V)  = 


K 

I 

n=l 


h  V 
n 


n 


(4.1) 
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In  this  case,  K  has  been  chosen  to  be  ten.  The  d.c.  charac¬ 
teristic  for  a  1N2940  tunnel  diode  has  been  obtained  graph¬ 
ically,  using  an  analogue  computer*,  and  is  shown  in  Fig.  4.1 
(a)  . 


If  ten  points  (V  ,1  ),  excluding  the  origin,  are 

n  n 

chosen  on  the  d.c.  characteristic,  ten  simultaneous  equations 

in  h  can  be  written, 
n 

i.e. 


hlVl  +  h2Vi 


+  h10Vl 


10 


(4.2) 


10 


hlV10 


+  h2V10 


+  h10V10 


10 


This  was  done  for  the  ten  points  indicated  in  Fig.  4.1(a), 
and  Eq.  4.2  were  then  solved  by  means  of  a  digital  computer. 
The  resulting  ten  coefficients  are  listed  in  Table  4.2 


Table 

Characteristic 

hx  =  +0.34887421xl0_1  A/V 
h2  =  -0.42478878  A/V2 

h3  =  +0.30580876x10  A/V3 
h4  =  -0.30233254xl02  A/V4 
hc  =  +0.26704247xl03  A/V5 

D 


4.2 

Coefficients 

hc  =  -0.14 33151 5xl04  A/V6 
b 

hy  =  +0.45227742xl04  A/V7 
hQ  =  -0.8312128 Ox 104  A/V8 

O 

h9  =  +0.82619199xl04  A/V9 
h10  =  -°-34370567xl°4  A/V10 


* 


Analogue  computer  simulation  is  discussed  in  All. 
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Figure  4 . 1 

Tunnel  Diode  Characteristics 


■ 
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Using  these  ten  coefficients,  the  approximate  function,  (4.1), 
was  evaluated  by  computer,  and  the  resulting  curve  is  shown 
in  Fig.  4.1(b).  This  is  seen  to  be  a  very  good  approximation 
to  the  experimental  curve.  Hence,  all  following  theory  shall 

be  based  on  the  above  values  of  h  . 

n 

4 . 2  Expansion  of  f (V)  About  the  Operating  Point 

It  will  later  be  necessary  to  have  an  expression, 

W (V) ,  such  that 

W (V)  =  f(V  +  V)  ;  (4.3) 

where  refers  to  an  operating  point  voltage.  This  is  il¬ 
lustrated  in  Fig.  4.2 


. 


23 


K 


Since  f (V)  =  £ 


n=l 


hnVn>  it  follows  that 


where 


K 


W(V)  =  l 


n=0 


g  V 


n 


'n 


f (n) (VQ) 
n  l 


W 


(4.4) 


(4.5) 


But 


f(n)(vQ)  = 


K  i! 

I 


h .  V, 


l-n 


i=n  (i-n) :  1  Q 


(4.6) 


Thus 


i! 


n 


K 

I  - 

i=n  n! (i-n) ! 


hiVQ 


i-n 


(4.7) 


Also,  by  definition. 


il 


n! (i-n) ! 


r  ■  ft  • 


(4.8) 


so  that  g  can  be  written  as 
n 


K  /  i' 


n 


l  U  hi 


v. 


i-n 


i=n 


(4.9) 


4 . 3  Basic  Circuit  Equations 

The  basic  oscillator  circuit  (steady-state  equi¬ 
valent)  is  shown  again  in  Fig.  4.3.  ZT  is  taken  to  be  the 
impedance  of  the  primary  tank  circuit,  and  Zfi  as  the  equi¬ 
valent  biasing  network  (Fig.  2.2). 


V) :  .i  2 


.t-i 


■ 


)  V  ~ 


irfj  3  _s  bi  e  ,ikuotio  *a*3  ^zmfxc  ■  ■  i  .  i  -  rtsl  imfc 

■  .  .  -  .  -•  •  ' 
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+  V 


D 


Figure  4.3 

Basic  Oscillator  Circuit 

Note  that  and  1^  denote  total  instantaneous 
D  D 

values  of  voltage  and  current,  and  V  is  a  fixed  d.c.  voltage. 

hi 

If  it  can  be  assumed  that  all  voltages  and  currents  are  La¬ 
place  Transformable,  -two  basic  equations  can  be  written  for 


the  circuit: 

Vs) 

=  VD(S)  +  ID(S) 

ZB(S)  +  zT  (S) 

(4.10) 

and 

ZD 

=  f<V 

(4.11) 

Now,  let 

VD 

-  VQ  +  V 

(4.12) 

and 

XD 

=  JQ  +  1 

(4.13) 

where  v  and  i  are  time-dependent. 

Since  the  circuit  is  known  to  oscillate  in  some 
manner,  it  can  be  assumed  that: 


' 
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and 


v  =  y  A  sin  (ncot  +  6  ) 
n=l  n  n 


=  Im 


00  j  (nwt+e  )~ 

V  A  -  n 


n=l 


n 


i  =  7  I  sm  nwt  +  <p  ) 

S  n  yn 

n=l 


=  Im 


00  j  (na)t+d>  ) 

T  1  -  n 

^  1 


n=l 


n 


(4.14) 


(4.15) 


For  d.c.,  ZB  =  RB  and  =  Rn ;  and  for  a.c.,  Zn  =  0  and 


B 


ZT  =  ZT  ( jnco)  . 

Upon  substitution  of  the  above  expressions  into 
(4.10),  and  letting  S  =  jnoo: 


VB  =  VQ  +  Im 


“  j  (no)t+6  ) 

l  A  e  n 

L  n 


n=l 


+  Iq(Rb  +  Rx) 


+  Im 


00  1  (ncot+d)  ) 

V  -r  n 

)  I  e 
L  -j  n 
n=l 


ZT(  jnco) 


(4.16) 


From  this,  it  is  seen  that: 


and 


VB  =  VQ  +  VR1  +  RB) 


J  0  l  <p 

•a  n  _  J  n  „  .  .  x 

A  e  =  -I  e  Zm(ina)) 

n  n  T  J 


(4.17) 


(4.18) 


If  the  impedance  of  the  tank  circuit  is  sufficient¬ 
ly  low  for  all  harmonics  of  the  frequency  of  oscillations,  it 


;  ; 
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can  be  assumed  that 


ZT(jnw)  =  0  ;  n  ^  1. 

Then,  from  (4.18),  it  follows  that 

j  6  T  j  4> 

A-^  x  =  -Iie  X  ZT(jo>)  ;  (4.19) 

and  =  0 ;  n  ^  1 .  Thus,  all  voltage  harmonics  are  neglected, 
and  the  analysis  is  confined  to  sinusoidal  or  near-sinusoidal 
oscillation . 

For  convenience,  0^  can  now  be  chosen  as  zero,  and 
1^  and  A^  can  be  denoted  simply  by  I  and  A.  Eg.  4.19  can 
thus  be  rewritten  as 

I  1 

-  - - tt  ;  (4.20) 

A  ZT  ( j  w)  e~*  ^ 

where  <J>  now  represents  the  phase  angle  between  the  fundamental 
of  current  and  the  fundamental  of  voltage. 

It  will  later  be  seen  that  in  neglecting  all  voltage 
harmonics ,  <f>  must  be  zero.  Hence,  from  (4.20),  ZT(jto)  must 
be  real,  or 

to  =  toQ  ,  and  ZT(jto)  =  Rq  (4.21) 

This  is  identical  to  the  result  derived  from  the  "Feedback 
Oscillator  Model". 

Equation  4.20  can  now  be  rewritten  as 
I  1 

R 

o 


A 


and 


(4.22) 


< 


■ 


■ 


27 


(4.17  can  be  rewritten  as 


I 


Q 


V 


B 


R 


B 


+  R. 


(4.23) 


Equations  4.22  and  4.23  are  the  two  equations  that 
must  be  simultaneously  satisfied  for  steady-state  sinusoidal 
oscillation . 


4 . 4  Describing  Function  and  Biasing  Equations 

Since  ip  =  I  +  i,  and  VD  =  Vq  +  v,  (4.11)  can  be 
rewritten  as 


(IQ  +  i)  =  f(V  +  v) 


(4.24) 


But,  from  (4.3),  f (Vq  +  v)  =  w(v)  ,  so  that 

K 

(I  +  i)  =  W(v)  =  l  gnvn 
w  n=0 


(4.25) 


Now,  since  v  =  A  sin  cat,  and  i  =  I  sin  cot  +  (Higher  Harmonics)  , 


it  follows  that: 


I  +  I  sin  cot  +  (Higher  Harmonics) 

K  _ 

r  n  .  n 

=  )  g  A  sin  cot  ; 

L  n  n 

n=0 


(4.26) 


from  which 


K 


Q 


component  of  £ 


n=0 


,n  .  n  j 

g  A  sin  cot,  and 


the  d . c . 


. 


03  %  (V)  W  *  (V  +  ..  %  ‘f  ,  ^ 
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K 


I  =  the  fundamental  component  of  £  g  Asinwt 

n=0  n 


The  following  formula  is  useful  for  the  expansion 


.  n 


of  sin  cot.  If  n  is  odd, 


•  n  I 

sm  cot 


(-1)  < n_ 1 > / 2  ^ 

-  I  <-Dr 

r=0 


,n-l 


sin  (n-2r)  cot  (4.27) 


Since  2r  is  even,  and  n  is  odd,  n-2r  must  be  odd.  Thus, 
sin11  cot  is  of  the  form: 


sm 


n 


tot 


a^(n)  sin  tot  +  a^(n)  sin  3cot  + 


...  +  a  (n)  sin 
n 

Equation  4.28  indicates  that  if  n  is  odd, 
fundamental  component  and  odd  harmonics. 

If  n  is  even,  n-1  will  be  odd, 
placed  by  n-1  in  (4.27).  Hence, 

n-1  (-D  (n“2) /2  V- 

sin  cot  =  - —2 -  l  (”!) 

2  r=0 


ncot .  (4.28) 

sin11  cot  has  only  a 


and  n  can  be  re¬ 


sin  (n-l-2r)  cot. 


Then , 

.  n 
sm 


.  n  .  n-1  , 

sm  cot  =  sm  cot  • 

(_l)(n_2)/2  SZZ 

cot  =  - - — o -  I 

2n'2  r=0 


sin  cot,  or 


sin  (n-l-2r )  cot  •  sin  cot 

(4.29) 


Since  n  is  now  even,  and  2r  is  even,  n-l-2r  must  be  odd,  so 
that  sin11  cot  is  of  the  form: 


■ 


-  ■  :  "  :  ' 


•  D  BCi  JftllfR  ybb  0  81  «  fc  ‘to  *  €»VB  '  >• 


r  :  >•  •  4  *  •  *  ,«-x 


■ 

ni6  ^Bfl 
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i*i  2 

sin  oot  =  aQ(n)  sin  wt  +  a2  ^  s^n  3cot  +  ... 


+  an_2  (R)  sin  wt  •  sin  (n-1)  oot 


(4.30) 


This  can  be  rewritten  as 


.  n  .  . 

sm  cot  =  a  (n) 

o 


(1  -  cos  2 oo t ) 


(cos  2cot-cos  4 cot) 


+  a2 (n) 


+  ...  +  a  0 (n) 
n-2 


[cos  (n-2 )  oot  -  cos  ncotQ 


(4.31) 


Thus,  if  n  is  even,  sin11  cot  has  only  a  d.c.  component  and  even 
harmonics . 


Referring  now  to  Eq's.  4.28  and  4.31  it  is  apparent 
that  I  depends  only  on  the  even  coefficients  g^,  and  the 
even  powers  of  An  sin11  cot;  and  that  I  depends  only  on  the 
odd  coefficients  gn  and  the  odd  powers  of  An  sinn  cot. 

First  consider  the  fundamental  component  of  cur¬ 
rent,  I,  as  given  by  (4.26).  The  general  term  in  the  expan¬ 
sion  for  (I  +  i)  is  gnAnsinn  cot.  The  fundamental  component 
due  to  this  term  (for  odd  values  of  n)  will  be  gnAna^(n), 

where  a-^(n)  is  the  coefficient  of  sin  cot  in  Eq.  4.27.  i.e. 

{_1)(n-l)/2  n-l 


a1(n) 


.n-l 


(-1)  11;  where  r  = 


Substituting 


n  =  2r  +  1  into  the  expression  for  a^(n); 

1  /  2r+l^ 

r 


a1(r) 


2r 


(4.32) 


. 


f 
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Also,  g^A31  sinn  ait  can  be  expressed  as  92r+lA^r+^S^n  a)t'  and  t^le 
fundamental  component  due  to  this  term  will  be 


, 2r+l  .  . 
g2r+i  A  ax(r)  = 


2r+l 


A 


2r+l 


2r 


(4.33) 


The  smallest  odd  value  of  n  is  one,  for  which 
2r+l  =  1,  or  r  =  0.  Also,  if  K  is  even,  the  largest  odd 
value  of  n  is  K-l ,  for  which  2r+l  =  K-l,  or  r  =  (K-2)/2. 
Thus,  the  total  fundamental  component  of  current  will  be 
given  by  the  equation 


I  = 


K-2 


r=0 


2r+l 


2r+l 


(4.34) 


Dividing  both  sides  by  A; 


(4.35) 


Equation  4.35  is  essentially  a  describing-function 
equation,  as  it  relates  the  amplitude  of  the  fundamental 
component  of  the  output  of  a  non-linear  device  to  the  am¬ 
plitude  of  an  applied  sinusoidal  input. 

K  /  i\  • 

Recalling  that  9=1  (  J  h^  V  ,  it  is 

n  i=n  '  '  y 

apparent  that  (4.35)  is  a  function  of  both  A  and  V^.  Further¬ 
more,  as  the  ratio  I/A  has  dimensions  of  conductance,  the 
summation  shall  be  defined  as 


G 

eq 


<A,Vq) 


(4.36) 


. 


•  1 


[A  vd  satiz  iJod  pnitiviO 


I 


' 
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Equation  4.22  can  now  be  written  as 


G  (A  V  )  =  - 

eq  '  Q 


(4.37) 


The  d.c.  component  of  current,  1^,  must  now  be 

found.  From  Eq.  4.31,  the  d.c.  component  of  the  general  term 

n  n  a  ^ 

g  A  sin11  ojt  (for  even  values  of  n)  will  be  q  An  — —  ; 

^n  ^n  2 

2 

where  aQ (n)  is  the  coefficient  of  sin'  wt  in  (4.29).  i.e. 


ao(n)  = 


,  ,  v  (n-2)  /2  /  i\  o 

(-1)  '  n-1  \  n-2 

- H=2 -  (-D  r  ''  where  r  = 


Substitut 


ing  n  =  2r+2  into  the  above  expression, 


ao(r)  = 


/2r+l\ 


:2r  V 


(4.38) 


Since  n  =  2r+2,  g  An  sin11  cot  can  be  written  as 

r  n 

2r+2  2r+2 

^2r+2  A  sin  wt,  and  t^le  ^*c*  component  due  to  this 


term  will  be 


g2r+2  A 


~  . o  a  (r) 
2r+2  o 


2r+2  2r+l  (2r+^\ 
=  g2r+2  A  (T}  V  r 


(4. 


The  smallest  even  value  of  n  in  the  expansion  of 
(1^  +  i)  is  zero.  This  corresponds  to  r  =  -1.  But  if  r  =  -1 

I'1] 

2r+l  =  -1,  and  the  expression  l  ^ J  is  meaningless.  By  in¬ 
spection  of  Eq.  4.25,.  it  is  seen,  however,  that  the  term  cor¬ 
responding  to  n  =  0  is  simply  gQ,  which  is  certainly  a  d.c. 
component.  i.e. 


K 

I 

i=0 


l 


hiVQ' 


K 

I 

i=l 


hiVQ' 


f(vQ) 


39) 
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The  next  smallest  even  value  of  n  is  two,  for  which  r  =  0o 
This  is  valid,  as  the  expression^  j  is  meaningful  0  Again 
assuming  K  to  be  even,  the  largest  even  value  of  n  will  be 
K,  for  which  r  =  (k-2)/2„  Thus, 


K-2 


g2r+2 


1 x  2r+l  _ 2r+2 


( 4  o  4  0 ) 


As  before,  gn  is  a  function  of  V^,  so  that  I  is 
a  function  of  both  VQ  and  A0  The  summation  in  ( 4  0  4 0 )  shall 
thus  be  defined  as 


XQ  (A'V 


(4.41) 


Equation  4e23  can  now  be  written  as 


VA'V  = 


RB  +  R1 


(4.42) 


It  is  shown  in  Sec0  5„6  that  a  solution  to  equation 
4023  (for  a  specific  value  of  A)  determines  and  I  «,  Hence 
this  equation  shall  be  referred  to  as  the  "Biasing  Equation" „ 
To  summarize  the  results  of  the  preceding  pages, 
for  steady-state  sinusoidal  oscillations,  the  following  two 
relationships  must  be  satisfied  simultaneously : 


G 


eq 


(A,V 


1_ 

R 

o 


<a'vq> 


+  R 


1 


( 4  o  43 ) 


' 
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Footnote 

It  was  stated  without  proof  in  Sec.  4.3  that  by 
neglecting  all  voltage  harmonics,  the  phase  angle  between 
voltage  and  the  fundamental  component  of  current  would  be 
zero.  This  is  true  because  in  the  expansion  of  sin11  cot,  for 
odd  values  of  n,  only  sine-terms  are  present.  A  phase-shift 
here  would  require  the  presence  of  cosine-terms  of  the  funda¬ 
mental  frequency. 

If  harmonic  distortion  is  permitted  in  the  voltage 
wave-form,  it  can  easily  be  shown  that  such  a  phase  shift 
can  occur.  For  the  sake  of  illustration,  it  will  be  suf¬ 
ficient  to  consider  only  the  second  harmonic,  and  the  second 
term  in  the  expansion  of  (1^  +  i) .  i.e. 

2  2 
=  g2  (A^  sin  cot  +  A2  sin  2cot) 

2  2 

=  g2  (Af  sin  cat  +  2A]_A2  s^-n  wt*sin  2ajt 
+  A^  sin^  2 cot)  . 

2  2 

Now  A^  sin  cot  contains  only  d.c.  and  second  harmonic,  and 
2  2 

A2  sin  2 cot  contains  only  d.c.  and  fourth  harmonic.  How¬ 
ever,  2A^A2  sin  cot  sin  2cat  =  A^A2  (cos  cat  -  cos  3cot)  .  Thus, 
the  term  g2A1A2  cos  cot  has  been  introduced,  which  is  a  cosine- 
term  of  the  fundamental  frequency.  This  must  produce  a 
phase-shift  between  current  fundamental  and  voltage  funda¬ 


mental  . 


r.  -«*i  '  ‘  '  '  ' 
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Referring  now  to  Eq.  4.20,  4>  will  no  longer  be 
zero,  so  that  ZT  ( j oi )  will  no  longer  be  real;  and  a  shift  in 
frequency  from  <*)  will  occur. 

It  must  be  concluded  then,  that  distortion  of  the 
voltage  wave  form,  in  general,  may  be  accompanied  by  a  de¬ 
parture  in  frequency  from  that  indicated  by  the  describing- 
function  method. 


ed  tepnol  on  Iliw  *  *0S,*  ..*>3  o*  won  paling  si  ^ 


V.  SOLUTION  TO  DESCRIBING  FUNCTION 


AND  BIASING  EQUATIONS 


5.1  Discussion  of  G  ( A , V  ) 
 eg  Q 


Because  of  the  complexity  of  Eq's.  4.43,  the  most 
convenient  means  of  finding  a  simultaneous  solution  is  to  de¬ 
velop  a  graphical  method.  The  initial  step  shall  be  a  rather 
detailed  discussion  of  the  describing-function,  . 

Using  the  ten  coefficients  listed  in  Table  4.2, 
the  describing-function  has  been  evaluated  by  means  of  a 
digital  computer.  Fig's.  5.1(a)  and  (b)  show  some  typical 
curves  corresponding  to  various  ranges  of  Vq. 

A  solution  to  Eq.  4.37  is  now  determined  by  the 

intersection  of  a  horizontal  line  at  -  1/R  with  the  curve 

o 

corresponding  to  the  value  of  chosen.  It  can  be  seen 

that  either  one,  two,  or  three  such  solutions  may  exist. 

The  function  G  (A,V~)  can  be  expressed  as 

eq  Q 


Geq(A'V 


gi  +  4  g3A  +  ( 


Higher  Powers 
of  A 


) 


(5.1) 


Thus , 


lim  G  (A,Vn) 
A-0  8q  Q 


9l  =  f(V) 


(5.2) 


Since  f'(V  )  is  defined  as  -  1/rd,  Eq.  4.37  now  becomes 


v  }o  l  fi  n  ii.  '*'©■  v  " 
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Figure  5.1  (a ) 
Describing “Function  Curves 
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Describing-Function  Curves 
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This  is  the  same  result  as  Eq.  2.9. 
From  (5.1) , 


9A 


3 

2 


.Higher  Powers^ 
g3  1  of  A  ' 


(5.3) 


and  hence , 


9G 


3SL 


9A 


-  0 


(5.4) 


A=0 


Thus,  all  describing  function  curves  meet  the  vertical  axis 
at  right  angles. 

As  a  curve  G  (A,Vq)  approaches  the  G^-axis,  it 
can  be  either  concave  upward  or  concave  downward.  This 
property  is  determined  by  the  sign  of 


92G 


eq 


3  A' 


But  from  (5.3) 


A=0 


so  that 


92G 


= 


+  (Higher  Powers,.  (5-5) 


9A 


„  g3  '  v  of  A 


92G 


eq 


9  A' 


=  -  g 


(5.6) 


A=0 


Thus,  if  g«(Vn)  >  0,  the  curve  is  concave  upward  at  A  =  0; 
and  if  g  (V  )  <  0,  the  curve  is  concave  downward  at  A  =  0. 


L  .  ■  ■  •>  J  d  ' 
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5.2  Radius  of  Curvature 


By  introducing  the  concept  of  the  radius  of  curva¬ 
ture  of  f (V)  at  the  point  V^,  a  relationship  between  the 
sign  of  g^  and  the  rate  of  change  of  the  radius  of  curvature 
(with  respect  to  arc-distance  along  the  curve)  can  be  derived. 
Note  that  in  the  Vl-plane ,  distance  will  have  mixed  units, 
(i.e.  S  =  /v  +  I  )  Thus,  all  physical  dimensions  shall  be 
neglected  in  this  section. 

The  radius  of  curvature,  p,  can  be  defined  as 


[l  +  f  (VQ) 2]  3/2 

f"<V 


(5.8) 


It  is  shown  in  Sec.  A1 . 1  that  the  rate  of  change  of  p,  with 
respect  to  arc-distance,  is  given  by 


dp  3f  '  (V)  f  '  '  (V)  2  -  [l  +  f  '  (V)  2]f  '  '  '  (V) 

dS  f''(V)2 


(5.9) 


The  basic  result  obtained  from  the  radius  of  curva¬ 


ture  concept  shall  be  simply  stated  here  as  a  theorem,  which 
is  proven  in  Sec.  Al . 1 . 

Theorem  5 . 1 

If  {X}  represents  the  complete  set  of  values  of 
VQ  for  which  f'(VQ)  <  0;  and  {R}  represents  a  con¬ 
tinuous  range  of  values  of  V^,  included  by  {X},  for 


which 


V=V, 


>_  0;  then  {R}  must  be  included  by  a 


■ 


■ 
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larger  continuous  range,  {R' }  for  which  g3 (vQ)  <  0. 
This  is  illustrated  in  Fig.  5.2. 


Graphical  Interpretation  of  Theorem  5.1 

It  should  be  pointed  out  that  if  f (V)  is  concave 
downward,  f''(V)  will  be  negative,  and  hence  p,  as  defined 
by  Eq.  5.8,  will  be  negative.  Now,  if  p  <  0,  dp/dS  <  0 
means  that  |p|  must  be  increasing.  This  is  the  case  in  the 
region  between  the  peak  point  and  the  inflection  point. 

As  the  relationship  between  the  sign  of  g^  and  the 
concavity  of  the  describing  function  curves  has  already  been 
discussed,  the  conclusion  to  be  drawn  at  this  point  is  that 
if  a  range  {R}  can  be  identified  (by  inspection  of  the  d.c. 
characteristic)  for  which  dp/dS  >_  0,  this  range  must  be 
characterized  by  describing  function  curves  which  are  con¬ 
cave  downward  at  the  G  -axis. 

eq 


£  -  a  :  - 
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: 

' 

■ 


41 


5.3  Constant-G  Chart 

_ 


If  ~  is  replaced  by  an  arbitrary  negative  con¬ 


stant,  r,  in  Eq.  4.37  the  resulting  equation. 


Geq(A'VQ: 


=  r 


(5.10) 


defines  a  locus  or  contour  in  the  AVQ-plane.  Hence,  choosing 

different  values  of  r  will  give  rise  to  a  family  of  such 

contours.  Such  a  family  of  contours  shall  be  referred  to 

as  a  Constant-G  Chart. 

eq 

Using  the  coefficients  listed  in  Table  4.2,  a 
digital  computer  was  used  to  numerically  find  a  solution  in 
A  to  Eq.  5.10,  Vq  and  r  being  specified.  A  sufficient  number 
of  values  of  Vq  and  r  were  chosen  as  to  obtain  reasonable 
accuracy  throughout  the  region  of  interest.  The  resulting 
chart  is  shown  in  Fig.  5.3. 

It  is  now  apparent  that  the  contour  Geg(A,VQ)  =  r 

represents  the  locus  of  all  possible  solutions  to  Eq.  4.37, 

if  r  is  chosen  to  be  -  i—  . 

o 

Note  that  the  Constant-G^  Chart  depends  in  no  way 
on  the  circuitry  external  to  the  tunnel  diode  itself,  but 
only  on  the  d.c.  characteristic  of  the  particular  diode  under 
consideration.  Such  a  chart,  if  obtained,  can  thus  be  used 
in  conjunction  with  any  desirable  type  of  network,  providing 
the  applied  voltage  is  sinusoidal,  or  nearly  so.  In  specifying 


n.  !o<  ria  , son;  H  .9  isXq-  /A  erfi  r.s  luomao  10  -tool  aeons 
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Figure  5.3 

Constant-G  Chart 
eq 
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a  tank  circuit  (i.e.  in  specifying  R  )  one  is  merely  confined 
to  a  particular  contour  on  the  chart. 


5.4  Properties  of  the  Constant-G  Chart 

 e<3 


It  was  shown  (by  Eq.  5.2)  that  lim  G  (A,V )  =  f ' (V  ) . 

A-0  eq  g  u 

Thus,  the  contour  G  (A,V_)  =  r  must  meet  the  V  -axis  at  the 

eq  Q  Q 

points  where  f'(V^)  =  r.  This  of  course  implies  that  the 

contour  G  (A,Vn)  =  0  must  meet  the  Vn-axis  at  V  =  V  and 
eq  v  U  u  * 

Vq  =  V^.  Furthermore,  it  implies  that  contours  for  which 

G  >0  must  meet  the  V.-axis  where  and  >  VTT;  and 

eq  Q  Q  P  Q  V 

that  contours  for  which  G  <0  must  meet  the  V^-axis  in  the 

eq  Q 

region  Vp  <  <  V^.  Also,  since  Geg(A,V^)  is  at  every  point 

a  single-valued  function,  no  two  contours  can  cross.  Thus, 

every  point  which  can  represent  a  solution  to  Eq.  4.37  must 

be  within  the  area  between  the  V^-axis  and  the  contour 

G  (A ,  =  0. 

eq v  '  Q 

If  the  function  G  (A,V,J  has  an  extremum  at  some 

eq  Q 

point, then  at  this  point 


9G 


=  0  , 


9G 


and 


-SSL  =  o  . 


9A 


3  V, 


3G 


Equation  5.4  shows  that 


®£L 


9A 


=  0.  Now,  under  partial 


A=0 


differentiation  with  respect  to  Vq,  A  is  considered  to  be  a 
constant,  so  that 


£  . 

. 


1  t- 
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9G 


33. 


3V. 


3V, 


A=0 


Geq(0'V- 


9G 


But  by  (5.2)  ,  Geg(0/VQ)  =  f  '  (V  )  .  Hence, 


=  0 


A=0 


implies  that  f  *  '  (V^)  =  0,  or  .  Thus,  the  point  (0,V.r) 

is  an  extremum  of  G  (7\,  V)  . 

eq  Q 

Also,  Geg(0,VI)  =  f'(Vj.)r  which  is  defined  as  -G^, 

by  Eq.  1.2.  This  is  the  minimum  value  of  f'(V^).  Thus,  it 

follows  that  no  point  can  exist  for  which  G  is  less  than 

eq 

-G^,  as  the  contour  through  such  a  point  would  have  to  meet 

the  V^-axis  where  f ' (V^)  were  less  than  -G^.  Hence,  the 

point  (0,V  )  must  be  an  absolute  minimum  of  G  (A,V^) . 

I  eq  Q 

It  can  now  be  concluded  that  all  contours  for 

which  Gg  <  0  must  enclose  the  point  ( 0 , ) ,  and  must  in 

turn  be  enclosed  by  the  contour  G  =0. 

J  eq 

It  was  pointed  out  earlier  that  for  certain  values 

of  Vq  and  Rq,  more  than  one  solution  to  the  equation 

G  (A,V^)  =  -  can  exist.  This  is  clearly  evident  from 

eq  Q  Rq  j 

the  Constant-G  Chart  as  well. 

eq 

It  is  shown  in  Sec.  Al . 2  that  the  contours  immedi¬ 


ately  surrounding  the  inflection  point  are  approximately 
elliptic;  with  centres  at  ( 0 , V-j. )  ,  major  axis  of 


'Ga +  r 

3g3 ( vx ) 


(5.11) 


avsrf  bi.iow  i  og  6  i  $••..!*«  •  ctnoc 
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and  minor  axis  of 


+  r 


(5.12) 


3g 3  (Vj.) 

The  major  axis  is  parallel  to  the  A-axis.  r  is  the  value  of 
G  on  the  contour,  and  is  of  course  negative; 


g3(VI}  = 


f ! ' ' (vI) 

3 ! 


and  Gd  =  -f ' (V^ 


For  the  diode  under  consideration,  and  g^V^) 

_  o 

have  been  evaluated,  and  are  approximately  7.36  x  10  A/V 

3 

and  0.545  A/V  respectively.  For  the  smallest  contour  shown 

-3 

m  Fig.  5.3,  r  =  -7.0  x  10  ,  so  that  the  minor  axis  should 

be 

/ (7.36  -  7.0)  x  10"3 

2  /  -  -  2  9.7  mv . 

/  3(0.545) 

This  can  be  easily  verified  by  inspection  of  Fig.  5.3. 


There  are  two  other  properties  of  the  Constant-G 


eg 


Chart  that  shall  be  stated  here  as  theorems,  and  are  proven 
in  Sec's.  A1 . 4  and  A1 . 5 . 

Theorem  5 . 2 

All  contours  G  (A,V_)  =  r  meet  the  V  -axis  at 

eq  Q  Q 

right  angles. 

Theorem  5 . 3 


If  a  contour  G  =  r  meets  the  V„-axis  at  a  value  of 

eq  Q 

Vq  within  a  continuous  range,  (Vq  },  for  which  all 

curves  G  (A,V^)  have  the  same  direction  of  con- 
eq  Q 

cavity  at  A  =  0,  then  the  direction  of  concavity 


' 
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of  the  contour  (at  A  =  0)  will  be: 

1)  the  opposite  to  that  of  the  curves  Ge^(A ,  V^) , 
if  Vq  >  Vj ;  and 

2)  the  same  as  that  of  the  curves  G  (A,\AJ  if 

eq  Q 


<  V 


I* 


This  is  illustrated  in  Fig.  5.4.  In  Fig.  5.4(a), 

V  and  V  belong  to  the  ranges  {V  }  and  {V  } 
t?i  u2  U2 

respectively,  as  shown  in  Fig.  5.4(b). 


G 

eq 


(A'V 


(a) 


Figure  5.4 

Graphical  Interpretation  of  Theorem  5.3 

Since  relationships  between  the  radius  of  curvature 
of  f (V) ,  g^ (Vq) ,  and  the  direction  of  concavity  of  the  de¬ 
scribing-function  curves  have  already  been  established,  it 
follows  that  these  relationships  can  be  extended  to  the 
Constant-G^  Chart.  Hence,  the  ranges  {R}  and  {R'},  as  shown 
in  Fig.  5.2,  would  be  characterized  by  contours  being  concave 


:  , 


' 
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upward  at  the  V^-axis. 

Figure  5.5  illustrates  the  basic  properties  of  the 
Constant-G^  Chart  that  have  been  discussed  in  this  section. 


Figure  5.5 

Basic  Properties  of  the  Constant-G  Chart 

eq 

One  minor  limitation  of  the  Constant-G  Chart  as 

eq 

derived  for  the  IN2940  tunnel  diode  should  be  pointed  out. 
Referring  to  Fig.  4.1,  it  is  seen  that  the  d.c.  characteris¬ 
tic  based  on  the  coefficients  listed  in  Table  4.2  is  known  to 
be  accurate  only  in  the  voltage  range  0  <_  V  <_  500  mv.  Thus, 
for  the  method  of  analysis  to  be  strictly  valid,  the  diode 
voltage  must  at  all  times  be  confined  to  this  range.  This 
means,  that  on  the  Constant-G^  Chart,  portions  of  contours 
within  the  regions  where  +  A  >  500  mv,  and  -  A  <  0  will 
probably  be  less  accurate  than  those  outside  these  regions. 
These  regions  are  indicated  in  Fig.  5.3. 


, 
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5 . 5  Discussion  of  the  Biasing  Function 

The  biasing  function,  I^(A,Vq) ,  based  on  the  ten 
coefficients  listed  in  Table  4.2  has  been  evaluated  using  a 
digital  computer,  and  the  resulting  curves  are  shown  in 
Fig.  5.6.  The  evaluation  was  done  only  within  the  range  of 
Vq  for  which  g^ (V  )  <  0 ,  as  this  is  the  only  range  for  which 

oscillations  may  occur. 

It  is  noticed  that  for  A  =  0,  the  biasing  function 
degenerates  to  the  d.c.  characteristic.  This  can  easily  be 
seen  from  Eq.  4.40.  i.e. 


lim  iq(a,v  )  =  gG  =  f(v  )  . 

A->0 

Also,  for  very  small  values  of  A,  1^  can  be  approxi¬ 
mated  by  the  first  two  terms  of  Eq.  4.40,  or 

Iq(A,Vq)  =  gQ  +  (5.13) 


Now,  it  is  known  that: 

1.  if  V  <  V  ,  g9  <  0  ;  f (V  )  is  concave  downward 

2.  if  V  =  V  ,  g9  =  0  ;  f (V  )  has  an  inflection  point 

3.  if  vq  >  Vj,  g2  >  0  ;  f (vq)  is  concave  upward. 

Hence,  from  (5.13),  it  follows  that  the  biasing  curves  cor¬ 
responding  to  very  small  values  of  A  will  be  beneath  f (Vq)  if 


V. 


if  V 


of  Fig. 


will  cross  f (VQ)  at  VQ  -  Vx ;  and  will  be  above  f (VQ) 
V  .  The  above  can  easily  be  verified  by  inspection 

5.6. 


'  V 


)«  no.,  4osI  nJ  rt  stri  ( ;  0  «  cp  ,TV  *  ~V  ix 
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Figure  5.6 
Biasing  Curves 
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506  Unique  Determination  of  A  and  V 


Referring  now  to  Eq0  4,42,  it  is  apparent  that  the 
quantity  (V^-V  ) / (R^+R^)  represents  a  straight  line  in  the 

i  a  slone  of  -  - - -  .  a  v 


V^Ig-plane;  with  a  slope  of  -  R  ■  ,  a  Vn-intercept  of  V 


B 


B 


and  an  I^-intercept  of  VB/(R^+RB)0  This  straight  line  shall 
be  referred  to  as  the  d.c.  load-line.  Thus,  the  solution  to 
Eq0  4,42,  for  a  particular  value  of  A,  is  determined  by  the 
intersection  of  the  d.c.  load-line  with  the  biasing  curve 
corresponding  to  that  value  of  A,  Fig,  4,7  illustrates  a 
set  of  biasing  curves,  along  with  a  number  of  possible  load- 
lines  o  Note  that  all  load-lines  shown  have  a  slope  of  - 

l/jil 

If  a  sufficient  number  of  points  of  intersection  are 
now  obtained  and  plotted  on  the  Constant-G^  Chart,  the  re¬ 
sulting  curve  represents  the  locus  of  all  possible  solutions 
to  Eq.  40420  Since  for  a  specific  load-line,  A  can  be  thought 
of  as  determining  V^,  the  functional  notation  =  V^(A)  can 
be  used  for  this  locus. 

If  it  is  now  recalled  that  the  contour  G _ (A,V^)  = 


eq 


Q‘ 


R 


is  the  locus  of  all  possible  solutions  to  this  equation, 


o 


it  is  apparent  that  a  point  of  intersection  between  this  con¬ 
tour  and  the  curve  Vq(A)  must  be  a  simultaneous  solution  to 
Eq'So  4,43,  Hence,  and  A  are  uniquely  determined. 

Figure  5,8  shows  the  curves,  V  (A)  corresponding 
to  the  load-lines  shown  in  Fig.  5.7 


. 


' 
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Figure  5.7 

Biasing  Curves  and  Load-Lines 
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Figure  5.8 
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Example  501 

As  an  example,  if  VB  were  chosen  as  225  mv,  and 

-  as  -4o0  x  10-3mho,  the  circuit  could  be  ex- 
o 

pected  to  oscillate  with  an  amplitude  of  111  mv0, 

and  a  of  152 „ 5  mv.  In  Fige  507,  this  corresponds 

to  the  third  load-line  from  the  left  and  in  Fig0  5C8 

it  corresponds  to  the  third  lowest  curve,  V^(A)0 

R  is  of  course  250fio 
o 

5 . 7  Stability  of  Solutions 

It  can  easily  be  seen  from  Fig0  508,  that  in  many 
cases,  there  is  more  than  one  solution  in  A  and  VQ  for  a 
specific  contour  and  a  specific  Vq(A)b  The  stability  of  the 
various  possible  solutions  will  now  be  analysed  and  it  will 
be  shown  that  both  "stable"  and  "unstable"  solutions  are  pos¬ 
sible  . 

A  "stable"  solution  shall  be  defined  as  a  solution 
for  which,  under  a  small  displacement  in  amplitude,  the 
system  returns  to  its  original  state . 

An  "unstable"  solution  shall  similarly  be  defined 
as  one  for  which,  under  such  a  displacement,  the  system  de¬ 
parts  from  its  original  state  and  tends  toward  some  other 
state  o 

Consider  the  linear  system  shown  in  Fig0  5090 
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Figure  5.9 

Comparative  Linear  System 


For  this  system 


zT(s) 


R  +  SL 


S2LC  +  SCR  +  1 


and 


Thus , 


Z  (S)  = 


I  (S) 
V(S) 


-  +  Z  (S) 
G 


(R  +  SL)G  +  S  LC  +  SCR  +  1 
G(S2LC  +  SCR  +  1) 


G(S  LC  +  SCR  +  1) 

S2LC  +  S (CR  +  GL)  +  (1  4-  RG) 


The  Characteristic  Equation  is 

S2LC  +  S (CR  +  GL)  +  (1  +  RG)  =  0  ,  and  its 


roots  are 


(CR  +  GL) 


S  =  - 


(1+RG)  (CR+GL) 


2  2 
4L  C 


(5.14) 


2LC  V  LC 

Now,  for  sinusoidal  oscillation,  the  roots  of  the 
Characteristic  Equation  must  be  imaginary,  or  (CR  +  GL)  =  0. 
From  this. 


CR 


G  =  - 


(5.15) 


R 


0  .  (.10  +  to)  to  9d  isaa  floiJ6Wpa  5  'U  *'"■ 
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Substituting  (5.15)  into  (5.14), 


S 


±j 


/l  -  (R2C/L) 
/  LC 


±jiu  (5.16) 
J  o 


Thus,  for  sinusoidal  oscillation,  the  describing-function  model 
of  the  non-linear  system  can  be  compared  to  the  linear  system 
having  a  negative  conductance,  G,  equal  to  Gecf  for  the  non¬ 
linear  system. 


Referring  again  to  the  linear  system,  if 


-1 


G  <  ,  -(CR+GL)  >  0 ,  or  the  roots  of  the  Characteristic 

o 

Equation  are  in  the  right-hand  side  of  the  S-plane.  This 
will  correspond  to  increasing  amplitude  of  oscillation.  Simi¬ 
larity,  if  G  >  ,  -(CR+GL)  <  0,  and  the  amplitude  will  de- 

JA 

O 

crease . 


Since  at  a  solution  for  sinusoidal  oscillation,  the 
describing  function  model  is  known  to  behave  in  a  similar  man¬ 
ner  to  the  linear  system,  it  can  be  assumed  that  the  same  will 
be  true  for  a  differential  displacement,  6A,  from  such  a  solu¬ 
tion.  In  the  AV^-plane,  the  displacement  6A,  corresponds  to 
displacement  along  the  curve  (A) ,  from  to  This  is 

shown  in  Fig.  5.10. 

The  displacement  6A  will  be  accompanied  by  an  in¬ 
cremental  change  in  G  ,  given  by 


dG 

6G  =  — ^  6 A 

eC*  dA 


(5.17) 


'  x  ■ 


ioao  l &bic  ua.  -  o.  to  it  .  o-:  £  it  i 


. 
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Figure  5.10 

Differential  Displacement  from  a  Solution 


Since  P^is  a  solution  for  sinusoidal  oscillation,  at  P^, 

Geq  =  ir  '  Hence'  at  P2 '  Geq  =  "  5“  +  31T2  6A’  If  ifc  is 
now  assumed  that  dG^/dA  is  positive  at  P^,  and  6A  is  taken 

1 


as  positive,  then  at  P^,  G  must  be  greater  than  - 

ec3  K0 

the  amplitude  should  decrease.  Taking  6A  as  negative,  and 


and 


dG  /dA  as  positive  would  result  in  G  being  less  than  —  , 
ec3  ec3  0 

and  the  amplitude  should  thus  increase.  Therefore,  it  is 
concluded  that  any  solution  for  which 


dG 


>  o 


(5.18) 


dA 


is  a  stable  solution.  A  similar  argument  shows  that  any 
solution  for  which 


dG 

— £2.  <  0  (5.19) 

dA 


is  an  unstable  solution. 


It  is  shown  in  Sec.  Al . 3  that  at  the  point  P^ 


:  ;  o  Xli  o  l  .  xurtxa  "toJ.  ro±3u  lo  a  i  rr.  aor.xB 
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dG  (r  -  r  )  cos(a-B) 

— ^  -  - - - —  — -  (5.20) 


dA  | AS |  cos  3 

In  Eq.  5.20,  r,  is  the  value  of  G  at  the  solution,  r~  is 

1  eq  2 

the  value  of  G  on  an  adjacent  contour,  such  that  r0  >  r,  ; 

eq  J  '  21 

and | AS |  is  the  perpendicular  distance  from  the  solution  to 
the  contour  G  =  Also,  3  is  the  acute  angle  between  the 

horizontal  axis  and  the  tangent  to  V^(A),  at  the  solution; 
and  a  is  the  angle  between  the  horizontal  axis  and  a  vector 
parallel  to  the  line  | AS | ,  in  the  direction  of  increasing  G  . 
This  is  shown  in  Fig.  5.11. 


A 


Figure  5.11 

Determination  of  Stability 

Now,  (r9-r, ) / | AS |  is  by  definition  positive,  and  since 

-90°  <_  3  <_  90°,  cos  3  must  be  positive  or  zero.  Thus,  the 

siqn  of  dG  /dA  will  be  the  same  as  the  sign  of  cos(a-3). 
r  eq 


■ 


, 
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Note  that  if  (a-3)  =  ±90°,  dG  /dA  =  0.  This  corresponds  to 

eq 

a  displacement  parallel  to  the  contour  G  =  r, .  Also  if 

eq  1 

6  =  ±90°,  dG  /dA  =  « .  This  corresponds  to  a  vertical  dis- 

eq 

placement,  so  that  dA  =  0 . 

Example  5 . 2 

Consider  now  an  example  where  there  are  known  to 

be  two  solutions.  Referring  to  Fig's.  5.7  and  5.8, 

this  could  correspond  to  VD  =  237.5  mv,  R, +RD  =  125ft, 

and  R  =  250ft.  There  is  seen  to  be  one  solution 
o 

for  A  =  116  mv  and  Vq  =  171  mv,  and  another  for 
A  =  58  mv  and  =  188.5  mv.  At  the  solution  for 
which  A  =  116  mv,  a  -  0,  and  6  a  0 .  Thus,  cos  (a-3) 

-  1 ,  so  that  this  will  be  a  stable  solution.  How¬ 
ever,  at  the  solution  for  which  A  =  58  mv,  a  >  90°, 
and  3  <  0°,  so  that  cos  (a-3)  <  0.  Thus,  this  solu¬ 
tion  will  be  unstable. 

This  means  that  if  the  oscillator  had  an  initial 

amplitude  of  less  than  58  mv,  oscillations  would 

die  out.  If,  however,  the  initial  amplitude  were 

greater  than  58  mv,  oscillations  would  build  up 

until  the  stable  solution  were  reached,  at  116  mv. 

It  can  also  be  seen  from  Fig's.  5.7  and  5.8,  that 

if  V  were  now  increased  to  250  mv,  there  would  be 
B 

no  solution;  and  that  for  some  value  of  Vg  between 
237.5  mv  and  250  mv,  the  stable  and  unstable 


. 
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solutions  would  converge  to  the  same  point.  This 
point  will  thus  determine  a  minimum  amplitude  for 
sustained  oscillation.  For  this  example,  the  mini¬ 
mum  amplitude  is  seen  to  be  approximately  105  mv. 
Note  that  oscillation  with  this  minimum  amplitude 

can  be  realized  only  by  gradually  increasing  V_ 

b 

from  a  smaller  value,  so  that  A  is  always  larger 
than  that  for  the  unstable  solution. 

A  similar  situation  would  occur  if  V_  were  decreased 

B 

to  somewhat  less  than  200  mv.  The  minimum  ampli¬ 
tude  in  this  case  would  be  about  65  mv. 

From  the  preceding  discussion,  it  is  apparent  that 
in  order  for  oscillations  to  build  up  from  zero,  the  solution 
occurring  for  the  lowest  value  of  A  must  be  stable  (if  more 
than  one  solution  exists) .  It  is  interesting  to  note  that  for 
VD  =  312.5  mv,  R, +R_  =  125ft,  and  -  —  =  -0.8  mmho ,  there  are 

B  _L  B  K 

O 

three  solutions.  The  first  and  the  third  are  stable,  and  the 
second  is  unstable. 

5 . 8  Relative  Stability 

Referring  once  more  to  the  linear  system  of  Fig.  5.9 
it  is  well-known  that  the  rate  of  growth  (or  decay)  of  ampli¬ 
tude  is  determined  by  the  magnitude  of  the  real  part  of  the 
root  S.  i.e.  the  magnitude  of  (CR+GL)/2LC.  If  this  is  also 
assumed  to  be  true  for  the  describing-function  model  under 


- 


* 
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the  virtual  displacement  6A,  it  then  follows  that  the  mag- 

•  (  T  2  —  T  i  \  cosffy.“R^ 

nitude  of  - - —  s  ^ —  will  determine  the  rate  with 

which  the  system  tends  toward  (or  away  from)  its  initial 

state .  Thus ,  the  quantity  cos  (a~e)  can  be  thought 

I  AS [  cos  3 

of  as  an  "Amplitude  Stability  Factor",  and  solutions  for 
which  this  quantity  is  relatively  large  can  be  termed 
"relatively  stable",, 

5 . 9  Further  Comments  on  Distortion 

Inherent  to  the  describing-function  method  is  the 
assumption  that  voltage  across  the  diode  is  purely  sinusoidal „ 

It  is  quite  obvious,  however,  that  the  current  waveform  can  have 
a  very  high  harmonic  content,  and  since  the  tank  circuit  does 
have  finite  impedance  at  any  frequency,  the  assumption  of 
sinusoidal  voltage  can  never  be  met  exactly.  It  was  also 
pointed  out  earlier,  that  for  large  degrees  of  voltage  distor¬ 
tion,  a  departure  in  the  frequency  of  oscillation  from  the 
resonant  frequency  of  the  tank  circuit  can  be  expected.  Al¬ 
though  a  detailed  analysis  of  distortion  would  be  too  diffi¬ 
cult  and  lengthy  to  be  practical, it  would  certainly  be  worth¬ 
while  to  determine  the  primary  factors  governing  distortion. 

Consider  once  more  the  "Feedback  Amplifier  Model" 
of  the  oscillator,  as  shown  in  Fig.  2.3(b).  It  will  be  re¬ 
called  that  an  equivalent  open-loop  gain,  Ay  =  -ZT/rd,  can 
be  defined.  This  is  valid  for  very  small  signals  only. 
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It  was  also  shown  in  the  discussion  of  the  tank 
circuit,  that  for  large  Q^,  the  impedance  to  the  n^  harmonic 
above  resonance  is  given  approximately  by 


ZT(jnu>o) 


Thus,  the  open-loop,  small-signal  gain,  at  a  frequency  of  n  , 


will  be  given  by 


Av(jnwo) 


n 


C1  (n  -l)rd 


From  this  it  is  apparent  that  if  (l/rd) /L^/C^  >  (n  -l)/n,  the 

open-loop  gain  for  the  nL  harmonic  will  be  greater  than  one 

Thus,  upon  closing  the  loop,  a  harmonic  initially  present  in 

the  voltage  wave-form  would  in  effect  be  amplified  and  fed 

back  to  the  input.  Hence,  it  must  be  concluded  that  if  the 

n1"  voltage  harmonic  is  to  be  negligible ,  (l/rd)  /L^/C^  must 

2 

be  much  less  than  (n  -l)/n. 


rd  ( i . e  . 
that 


If  the  second  harmonic,  and  the  minimum  value  of 
1/G-,)  are  considered,  the  resulting  relationship  is 


B..I 


^ /  c 

U1 


(5.21) 


It  is  to  be  expected  that  voltage  distortion  will 
be  decreased  as  amplitude  is  decreased,  as  the  harmonic  content 
of  the  current  wave-form  will  be  decreased.  Thus,  for  small 
values  of  amplitude,  near-linear  oscillation  may  be  realized 
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even  though  (5.21)  is  not  satisfied. 

Referring  now  to  Fig.  3„lp  it  has  been  shown  that 
if  /l^Tc^  is  increased: 

1»  Q-^  will  be  increased; 

2C  the  impedance  to  all  harmonics  will  be  increased;  and 
3.  Rq  will  be  increased. 

In  general,  increasing  Rq  will  result  in  an  increase  in 
amplitude,  and  thus  an  increase  in  the  harmonic  content  of 
the  current  wave-form.  Now  since  the  impedance  to  all  har¬ 
monics  will  be  simultaneously  increased,  it  follows  that  an 
increase  in  voltage  distortion  may  occur,  in  spite  of  the 
fact  that  will  have  been  increased. 

It  is  of  interest  to  note  that  for  a  tank  circuit 
having  R^  =  0,  both  and  Rq  would  be  infinite.  Thus,  vary¬ 
ing  /L^/C1  would  do  nothing  more  than  produce  different  de¬ 
grees  of  voltage  distortion  and  frequency  departure.  Since 
an  infinite-Q  tank  circuit  can  easily  be  simulated  on  an 
analogue  computer,  this  suggests  an  interesting  means  of  ex¬ 
perimentally  verifying  the  previous  conclusions, (See  SecoA203) 

5.10  Effect  of  Self-Oscillation 

It  was  stated  in  Sec.  1.3  that  normally  a  tunnel 

diode,  when  used  in  a  circuit,  will  have  associated  with  it 

parasitic  circuit  elements,  R  ,  L  ,  and  C  ;  and  an  equi- 

s  s  s 

valent  circuit  for  a  real  tunnel  diode  was  suggested.  This 
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Equivalent  Circuit  for  a  Real  Tunnel  Diode 

If  at  some  instant  a  d.c.  voltage  appears  across 
the  terminals  of  the  real  diode,  as  shown  in  Fig.  5.12,  then 
for  high  frequencies,  the  two  terminals  are  essentially  con¬ 
nected  together.  Thus,  the  ideal  diode  is  connected  across  a 
parallel  resonant  circuit.  It  then  follows  that  if  the 

resonant  resistance,  L  /R  C  ,  is  larger  than  1/Gd  for  the 

s  s  s 

ideal  diode,  self-oscillation  will  occur  whenever  VD  is  such 

that  the  ideal  diode  is  appropriately  biased.  If  R  is  quite 

s 

small,  it  can  be  neglected  for  d.c.,  so  that  the  d.c.  voltage 
across  the  ideal  diode,  V^,  will  be  the  d.c.  voltage,  . 

If  it  is  assumed  that  self-oscillation  of  amplitude 
does  occur  for  some  value  of  VD  (say  VD  )  the  d.c.  current 
through  the  real  diode,  I  ,  can  be  found  directly  from  the 
Biasing  Curves,  by  taking  as  VD  and  A  as  A^ . 

If  some  value  is  assumed  for  the  parasitic  resonant 
resistance,  Ls/Rscs'  the  contour  Geg  =  -RgCs/Ls  will  deter¬ 
mine  the  amplitude  of  self-oscillation  for  any  d.c.  voltage, 


- 

■ 


s  t'  :  5  .V  Vd  *  '  -Wl>  izsi  3 


)  \  •  ,  ■>  '  •  - 
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VQ'  across  the  ideal  diode.  Since  it  has  been  assumed  that 

Rg  can  be  neglected  for  d.c.,  the  load-line  for  the  ideal 

diode  will  be  vertical,  so  that  the  curve  V  (A) ,  on  the  Con- 

stant-G  Chart,  will  be  the  straight  line  V  =  V  .  This  is 

Q  D 

shown  in  Fig.  5.13. 


4 


Figure  5.13 

Derivation  of  the  Modified  D.C.  Characteristic 


For  a  contour  such  as  that  shown  in  Fig.  5.13  there 

can  again  be  two  solutions.  As  before,  the  one  for  the  smaller 

value  of  A  will  be  unstable,  and  the  other  one  will  be  stable. 

Thus,  if  V  were  increased  from  V  to  V  ,  A  would  increase 
^  U1  ^4 

from  zero  to  a  maximum,  and  then  decrease  until  V  were 

D4 

reached.  At  this  point,  oscillation  would  of  course  cease 

completely.  If  VD  were  now  decreased,  oscillation  would  not 

re-occur  until  V  were  reached,  and  the  amplitude  would  then 

°3 

vary  continuously  until  V  were  again  reached. 

1 

It  is  now  apparent  that  if  a  number  of  points  on 

=  -L  C  /R  are  plotted  on  the  Biasing  Curves 
eg  s  s'  s  c 


the  contour  G 


* 


rtxg  . sboib  Isebi  ad*  as 


fV  ad 


c  '  c  •  •  -  -V 


■ 
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a  Modified  D.C.  Characteristic"  will  result.  This  modified 
characteristic  is  analagous  to  that  defined  in  the  Dual  In¬ 
put  Describing-Function  method*,  with  the  exception  that  the 
self -oscillation  is  present  only  for  a  certain  range  of 
values  of  V^,  and  within  this  range  varies  in  amplitude  be¬ 
tween  zero  and  some  maximum  value.  Thus,  the  d.c.  characteris¬ 
tic  will  be  modified  only  in  the  range  where  self-oscillation 
occurs,  and  the  degree  of  modification  will  depend  on  the  am¬ 
plitude  of  oscillation. 

Referring  again  to  Fig.  5.13,  for  increasing  VQ,  the 

characteristic  will  be  modified  everywhere  in  the  range 

VD  <  VD  <  VD  ;  and  for  decreasing  VD,  it  will  be  modified  in 

the  range  <  v_  .  Thus,  the  modified  d.c.  character- 

D1  D  D3 

istic  will  exhibit  hysteresis. 

Assuming  a  parasitic  resonant  conductance  of 

_3 

3.0  x  10  mho,  a  modified  d.c.  characteristic  has  been  de¬ 
rived,  and  is  shown  in  Fig.  5.14.  This  form  of  characteris¬ 
tic  is  often  seen  when  a  tunnel  diode  characteristic  is  dis¬ 
played  on  a  conventional  curve  tracer  (see  Ref.  3,  page  67) . 

Consider  now  the  real  tunnel  diode  used  in  con¬ 
junction  with  an  external  tank  circuit  and  biasing  network  as 
shown  in  Fig.  5.15 


* 


See  Ref.  1,  page  408. 
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6 
i— i 

Figure  5.14 

Modified  D.C.  Characteristic 


AUJ  ' 
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R 


1 


Figure  5.15 

An  Oscillator  Using  a  Real  Tunnel  Diode 


If  it  is  assumed  that  the  resonant  frequency  of  the  parasitic 
tank  circuit  is  much  higher  than  that  of  the  external  tank 
circuit,  then  the  terminals  of  the  real  diode  will  again  be 
essentially  connected  together  for  voltages  at  the  frequency 
of  self-oscillation .  Furthermore,  during  any  one  cycle  of 
the  self-oscillation  (if  it  does  occur)  VD  can  be  considered 
as  relatively  constant,  so  that  the  modified  d.c.  characteris¬ 
tic  can  be  used  to  find  the  corresponding  current  through  the 
real  diode.  Hence,  the  circuit  of  Fig.  5.15  can  be  analysed 
by  the  methods  discussed  in  this  thesis,  except  that  the 
modified  d.c.  characteristic  must  be  used  in  place  of  the 
curve  I  =  f (V) .  This  of  course  means  that  if  the  describing- 
function  method  is  to  be  used,  a  modified  Constant-G  Chart 

tiq 

and  modified  Biasing  Curves  would  have  to  be  derived  on  the 
basis  of  the  modified  d.c.  characteristic.  A  much  more 


a  3  ,!J  jn  nli  oq. iMO,  «U  ^  IS  :  ■'  •«  *  «*>  >** 
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practical  approach  is  to  avoid,  if  possible,  conditions  lead- 

% 

ing  to  self-oscillation. 

In  conclusion  ,  the  describing-function  method,  as 
outlined  in  this  thesis,  can  be  expected  to  yield  reasonably 
accurate  results  only  if  conditions  are  such  that  voltage  is 
nearly  sinusoidal,  and  if  the  physical  layout  of  the  circuitry 
is  such  that  self -oscillation  does  not  occur. 
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VI.  Effect  of  Second  Tuned  Circuit 

It  has  been  established  that  the  circuit  shown  in 
Fig.  1.4(a)  can  be  made  to  oscillate  with  frequency  w  ,  such 
that  Im  { Z T  ( j  03^ )  }  =  0,  providing  a  stable  solution  to  Eq  '  s  . 

4.43  can  be  found.  Furthermore,  aside  from  the  Biasing  Curves 
and  Constant-G^  Chart,  which  pertain  only  to  a  particular 
tunnel  diode,  the  determination  of  a  solution  to  Eq '  s .  4.43  re 
quires  only  the  specification  of  the  d.c.  load-line,  and  the 
resonant  resistance  of  the  tank  circuit.  Since  mutually  coup¬ 
ling  a  second  tank  circuit  cannot  alter  the  d.c.  load-line,  it 
follows  that  analysis  of  this  system  will  require  only  addi¬ 
tional  knowledge  concerning  its  resonant  states. 


6 . 1  Resonant  Frequencies 


In  the  following  analysis,  ZT(S)  shall  refer  to  the 
impedance  of  the  mutually  coupled  tank  circuits,  as  shown  in 
Fig .  6.1;  or 


V(S) 

>  /  *>  »<Mpr 

ZT(S)  =  -  .  (6.1) 

I  (S) 


Figure  6.1 

Mutually  Coupled  Tank  Circuits 
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where 


It  is  not  difficult  to  show  that 


ZT(S)  = 


_  3  2 

a_S  +  a0S  4-  a, S  +  a 
3 _ 2. _ 1 _ o _ 

b.S4  +  b-.S3  +  bnS2  +  b,  S  +  b 

4  3  2  1  0 


(6.2) 


o 


=  R. 


L1  +  R1R2C2 
R1L2C2  +  L1R2C2 
L1L2C2(1  '  r2) 


o 


(6.3) 


t  R2^2 

L-^C^  +  3L» 2 C 2  t  ^i^2<31<32 

R1R2<31^'2  R1R2(31<32 
L1L2C1C2(1  -  r2)  ' 


and 


K‘ 


M‘ 


L1L2 


.  Note  that  0  <  K  <  1 


(6.4) 


Letting  S  =  jw, 


ZT(j<*>) 


2  .  3 

(aQ  -  )  +  j  (a^oo  ~  a3w  ) 

T 


4  3 

(bQ  -  b2w  +  b4co  )  +  j  (b^oj  -  b^co  ) 


(6.5) 


For  resonant  frequencies,  Im{ZT(jo))}  =  0,  which  results  in 
the  equation: 

6  4  2 

a)  a^b,,  +  a)  (b2^2  ~  a3^2  "  a1^4^  +  00  +  aib. 


■3  4 


’3  0  12 


-  bxa2  -  b3a0)  +  (aQbx  -  a±b0)  = 


0 


(6.6) 


<■><  £•  ~  r  :  >  '  o'  .  1 
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•  •  2 
Equation  6.6  is  a  cubic  in  oo  ,  and  cannot 

solved  without  loss  of  generality.  However,  it  may 

either  one,  two,  or  three  real  solutions  in  oo^. 

Note  that  if  (a  b.,  -  a,b  )  =  0,  a  trivial 

o  1  1  o' 

2 

to  (6.6)  is  oo  =0.  Substituting  the  values  of  aQ , 
and  bQ  into  the  above,  the  result  is  that  Q  ^  =  1. 


readily  be 
have 


solution 


i .  e . 


aobi  =  Ri<Rici  +  R2C2)  =  Rl2ci  +  R1R2C2  ;  and 

ajbo  =  +  R..R2C2  ,  so  that 


Since  for  a  single  tank  circuit. 


to 


o 


Llcl 


(1  -  — j)  ,  the 


above  corresponds  to  the  case  where  the  primary  circuit  is  not 


resonant . 


Because  of  the  inconvenience  that  would  be  involved 
in  analytically  solving  Eq.  6.6,  a  numerical  means  of  deter¬ 
mining  the  resonant  frequencies  will  now  be  developed.  This 
method  will  also  permit  a  certain  amount  of  optimization  of 
the  resonant  states  if  it  be  desired. 

Consider  the  circuit  shown  in  Fig.  6.2(a),  where 

V  (S) 

Z,  (S)  =  — - 

1  IX(S) 


(6.7) 


' 


- 
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I1(S)  I1(S) 


Figure  6.2 
Equivalent  Circuits 


Analysis  shows  that 

{S2L  C9(l-K2)  +  SR9C  +  1} 

Z  (S)  =  SL  - - — - 

(s  l2c2  +  sr2c2  +  1) 

and 

{1  -  (jo2L2C2  (1-K2)  +  jwR2C2> 

Z,  (jw)  =  jwL,  - 9 - 

(1  -  a)  L2C2)  +  jwR2C2 


Equation  6 . 9  can  be  written  as 


where 


(a1  +  j  3) 

Z  ( j  a) )  =  jo)L,  -  , 

(a2  +  je) 


a,  =  1  -  oj2L2C2(1-K2) 

2 

a2  =  1-a)  L2C2  ,  and 

3  =  o)R2C2 


(6.8) 


(6.9) 


(6.10) 


(6.11) 


3 


. 


73 


Rationalizing  (6.10)  results  in  the  following  expression: 


coL^3(a^  -  a0)  ^ala2  +  ^  ^ 

_  +  j  co  L  - 

(a2Z  +  3Z)  1 


2  2 
(a2Z  +  3  ) 


(6.12) 


Substituting  (6.11)  into  (6.12); 


Z1(ja)) 


4  2  2 

“  L1L2C2  R2K 


j  ooL. 


+ 


2  2  2  2  2 
(1  -  03  L?C7)  +  0)  R2  C2 

2-  -  —  2-  -  ”  -2”  ■  2-  2-  2] 


_(l-03  L2C2)  {1-03  L2C0  (1-K  )}  +  0)  R0  C7 


7  2  2  2  2 

(1  -  03  L2C2  )  Z  +  03  R2  C2 


(6.13) 


If  Im  {Z,  ( j 03 )  }  is  defined  as  jo3L  ,  and  Re{Zn(jo3)} 

1  J  J  03  1  J 


as  R  ,  (6.13)  can  be  written  as 


03 


Z  t  ( j  to )  =  R  +  jo3L 

1  J  03  J  03 


(6.14) 


Thus,  it  can  be  concluded  that  the  two  circuits  in  Fig.  6.2 
will  have  identical  impedances. 

Normalized  frequency,  u,  can  now  be  defined  such  that 


2  2 

u  =03  I*2C2  ; 


(6.15) 


and  Q0  such  that 
2 


Q. 


R22c2 


(6.16) 


Note  that  u,  as  defined  above,  has  a  different  meaning  than 


in  Eq.  3.5. 

Also,  let  R 


R 


03 


n 


and 


(6.17) 


R. 


03 


n 


(6.18) 


. 
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Substituting  (6.15)  and  (6.16)  into  (6.14),  the  results  are: 

2  -i 


u 


0) 


n 


IJl-u2)  {1-u2  (1-K2)  }  +  Q  2  _ 
. .  _  2 

9  7  7  7 

{  (l-u)z  +  uVQ2  > 


(6.19) 


and 


R 


R 


0) 


4  2 
u  K 


L,  R~ 


n 


R. 


/  / 1  2^2  2  ,  2  x 
{ (1-u  )  +  u  /Q0  } 


(6.20) 


The  quantity  shall  now  be  defined  as 


L2C2 

L1C1 


(6.21) 


This  is  approximately  the  ratio  of  the  resonant  frequency  of 
the  primary  tank  circuit  to  that  of  the  secondary.  Also 
shall  be  defined  as 


(6.22) 


Consider  now  the  ratio  Qe/C>2'  as  given  by  (6.22), 
(6.21)  and  (6.16)  . 


Now,  by  (6 . 23)  ,  (6 . 20)  can  be  written  as 


R 

n 


u4k2 (Qe/Q2> 

2  2  2  2 

(l-u)Z  +  uZ/Q2Z 


(6.24) 


, 
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If  all  previous  definitions  are  kept  in  mind,  the 
circuit  of  Fig.  6.1  can  now  be  represented  by  the  one  shown 
in  Fig .  6.3. 


Figure  6.3 

Equivalent  Circuit  for  Mutually  Coupled  Tank  Circuits 


For  this  circuit, 


YT(jw)  = 


ZT(jo>) 


=  jooC,  + 


(R,  +  R  )  +  j  ooL 

1  00  J  00 


(R,  +  R)  —  j  ooL 

^  t  1  00 
jooC-.  t  o  j 

1  (Rn  +  R  )  Z  +  oo  L 

1  00  00 


.  (6.25) 


For  resonance,  Im{YT(joo)}  =  0,  or 


oo 


2  2  2 

(R,  +  R  +  ooZL 

1  00  00 


(6.26) 


But 


L  =  L, L  ,  and  R  =  R, R  ,  so  that  (6.26)  can  be  written 
oo  1  n'  oo  In 


as 


' 


30  ,0  =  {<«t)TX»«I  ,9t.aenoaai  io« 


76 


L,L 
1  n 


1  2  M  _  .  2  ,  2  2  2 
Rn  (1  +  R  )  +co  Ln  L 
1  n  In 


L,L 
1  n 


r12U1  +  Rn)2  +  a.2Ln2  (L1/R1)  2} 


This  can  be  written  as 


r12c1 


2  2  2  2 
(1  +  R  )  +  co  L  Z  (L,/R,  ) 

n  nil 


n 


But 


r12c1 


Q- 


,  so  that 


(1  +  R  )  2  +  co2L  2  (L,/R,  )  2 
n  nil 


(6.27) 


n 


2  2 

By  Eq.  6.12,  the  term  oo  (L^/R^  can  be  expressed 


as  (u2/L2C2j  (L1/R1)2. 


2  2 


Consider  now  the  ratio  /F ^  .  Since  = 


Ll/Rl2Cl,  and  Ff2  =  L2C2^L1C1  ' 


L1  L1C1 


R1  C1  L2C2 


L1  2 

-  (— )  . 


L2<22  R1 


Thus , 


2 /Li\ 2 
00  ( - ) 

R, 


=  u 


2  Qi 


.  But  from  (6.22) 


0  2/f  2  =  o  *  ,  so  that  (6.27)  can  now  be  written  as 

1  '  f  e 
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2  2  2  2 
( 1  +  R  ;  +  u  Q  ZL  * 

n  e  n 


(6.28) 


Substituting  (6.19)  and  (6.24)  into  (6.28),  the  resulting 
equation  is 


{  (1-u2) 2  +  u 


Q- 


Q. 


42  222  2  2  2  u22 

+uqKz  ~^}Z+uQ  Z{  (1-u^)  {1-u^  (1-IT)  }  + 

2  Q2  e  Qo" 


2>2  u2  2 .  2/n  ^2*  ,  u“ 

{(1-u  )  +  — ~}{(l-u  )  {1-u  ( 1-K  )}  +  — 


(6.29) 


The  solution  to  Eq.  6.29  now  determines  the  resonant  frequen¬ 
cies  to  the  circuit  shown  in  Fig.  6.1. 

The  right-hand-side  of  Eq.  6.29  shall  be  defined  as 

2  2 

F(u  ,K  ) ,  so  that  (6.29)  can  be  written  as 

Q12  =  F(u2,K2)  (6.30) 

It  should  now  be  pointed  out  that  there  are  only 

2  2 

two  independent  parameters  in  F(u  ,  K  )  .  i.e.  anc^  Qe* 

If  Q9  and  Q  are  specified,  and  the  function  is  evaluated  and 

m  0 

2 

plotted,  u  ,  (for  resonance)  will  be  determined  by  the  mter- 

2  . 

sections  of  the  line  at  with  the  function.  Thus,  in 

varying  ,  the  solutions  can  be  shifted  one  way  or  the 
other  as  desired.  It  must  be  kept  in  mind,  however,  that  in 
varying  Q1,  Qe  must  remain  constant,  and  since  Q&  =  Q1/Fff 
it  follows  that  Ff  must  be  varied  correspondingly. 


■ 


f,i£  *  (.|SM)VinVl)JI^  - 


.I.i  .pi*  n X  nwoAK  31UOT.1  o  srf*  od  esi- 

-  X 

L  1 _  8£  nsdiiiw  9d  neo  (6£.d>  derfd  oa  vi  ••*<»  ' 


■ 
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2  2 

6.2  Discussion  of  F (u  ,K  ) 


Since  Eq.  6.30  determines  the  resonant  frequencies 
2  2 

for  any  value  of  K  ,  it  must  also  apply  for  K  =0,  which  cor- 

2 

responds  to  a  single  tuned  circuit.  Letting  K  =0, 


,  2  ,  2  2 
F  (u  ,  0 )  =  1  4-  Q  U 


(6.31) 


Thus ,  from  (6.30)  , 


2  2  2 
1  +  Q  u  =  Q, 
e  1 


or 


u 


Qi  -  1 


Q. 


But  Q 


2  2  ve 

(Qz-  1)fJ  e2 

so  that  u  =  - 2 -  =  Ff  (1  ” 


(6.32) 


2  2  2 
Now  u  =oo  L2C2  ,  and  F^ 


expressed  as 


L2C2 

L1C1 


Thus,  (6.32)  can  be 


w  ^2^2 


L2C2 

L1C1 


(1 - j)  ,  or 

Q  i 


00  = 


/Lici 


/  a. 


2  * 


This  is  the 


resonant  frequency  for  a  single  tank  circuit. 

In  the  following  discussion,  only  positive  values 
o 

of  u  are  considered. 

Referring  again  to  Eq.  6.31,  it  is  seen  that  for 
Y?  =  0,  F(u^,K^)  is  a  straight  line,  with  a  vertical  intercept 


. 


.  1  :■*  ■  ■  ■  • ■  , 
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( 


2 

of  one,  and  a  slope  of  .  Also,  it  can  be  shown  that 
lim  F(u2,K2)  =  1,  and  lim  F(u2,K2)  =  °°. 

u^  +  0  u2->oo 

Further  inspection  of  (6.29)  indicates  that  the 
2  2 

numerator  of  F(u  ,K  )  must  always  be  greater  than  zero,  be¬ 
cause  the  first  term  must  always  be  positive,  and  the  second 
term  must  always  be  positive  or  zero.  Furthermore,  the  first 

term  in  the  denominator  must  always  be  positive.  Hence,  the 
2  2 

sign  of  F(u  ,K  )  will  be  the  same  as  that  of  the  second  term 

2  2 

in  the  denominator,  and  F(u  ,K  )  will  have  vertical  as^ymptotes 
wherever  this  term  is  zero?  or 


(1 


2. 

u  ) 


2  2 
u  (1  -  K  ; 


0  (6.33) 


Equation  6.33  can  be  rewritten  as 


u4(l  -  K2)  -  u2(2  -  K2  - T)  +1=0 


(6.34) 


The  roots  of  this  are 


u 


(2  -  K2 - fr)  ±  /(2-K2-  1/Q,2) 2-4(l-K2) 

_ Q 2Z _ I _ 

2(1  -  K2) 


(6.35) 


Consider  first  the  case 
as^ymptote,  or  one  root  to  (6.34). 
4(1  -  K2) ;  which  can  be  written  as 


for  which  there  is  only  one 


+ 


0 


' 
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From  this, 


K‘ 


Q- 


T  V 


«2* 


+  4  (- 


Q2  $2 


1  2 
— _  +  _ 


Qo  Q 


or 


K‘ 


Qo  Q' 


Substituting  (6.36)  into  (6.35) 

2 

Q 


2  -  —  +  — 


u 


Q27  q22 


2(1  -  —  +  -ir) 


(6.36) 


,  which  simplifies 


to 


u 


°2  -  1 


(6.37) 


The  final  location  of  the  two  as^ymptotes ,  as 
I^+l,  can  easily  be  found.  Equation  6.35  can  be  written  as 


u 


N1(K2) 

D(K2) 


or  u 


n2(k  ) 

D(K2) 


;  where : 


(1)  lim  N,  (K  ) 

k2+i  1 


11  1 
=  (1 - y)  +  (1 - y)  =  2(1 - y) 

Qo  Qo  Qo 


(2)  lim  N0(K  ) 
K2^1  2 


=  (1  -  — y)  -  (1  -  — y)  =  0;  and 

Qo  Qo 


lim  D  ( K  )  =  0  . 

K2-*i 


(3) 


. 
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Thus , 


lim 


N, (K2) 


K2+l  D(K2) 


(6.38) 


By  the  application  of  1' Hospital's  Rule, 


1  im  - yr— =  1  jjfl  “ 

K2-l  D(K2)  K2^l  D' (K2) 


k2 (K2) 


n2' (IT) 


N '  (K“)  =  - T 

dK 


(2  -  K‘ 


- ~)  -  J ( 2  - 


2  2  2  2 
K  -  1/Q2)  -4(1-0 


=  -1  - 


-2(2  -  K2  -  1/Q22)  +  4 
2/(2-K2-l/Q22)  2  -  4  ( 1— K2 ) 


Similarly 


D'(K2)  = 


dK2  L. 


2 ( 1— K  / 


=  -2 


Thus , 


N  ' (K2)  1  1  2(2  K“q2)+4 

_2 _  =  _  +  _ 22 _ 

D' (K2)  2 


4  / (2-K2-l/Q0^) ^  -  4(1-Kfc) 


2 .  2 


Letting  =  1,  and  simplifying, 


N2'(K  )  2 

lim  - « —  =  u 


Q. 


K2+l  D' (k2) 


V  - 1 


(6.39) 


The  preceding  discussion  can  now  be  summarized  by 


the  following  properties: 
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(1) 

lim 

K2+0 

2  2 
F  (u  ,  K  ) 

(2) 

lim 

u2+0 

2  2 

F (u  ,K  ) 

(3) 

lim 

■Q  2  — oo 

2  2 

F  (u  ,  K  ) 

(4) 

2 

if  K 

2 

1  +  Qe2"2  ; 


=  1  ; 


1 


2  „2 


Q  o  Q  • 


2  f  F(u  ,K  )  is  positive  and  continuous; 


(5)  if  K‘ 


2 

Q 


2  ^2 


2  2 

F(u  , K  )  is  positive 


Q. 


2  ^2 

and  has  one  vertical  as/ymptote  at  u  = 


q2  -  1 


2  2  1  2  2 

(6)  if  K  >  —  -  — j  ;  F (u  , K  )  has  two  vertical  as^ymptotes 

°2  °2 

(given  by  Eq.  6.35)  between  which  it  is  negative. 
Everywhere  else  it  is  positive  ;  and 


(7)  as  K  >1,  these  as^ymptotes  approach  +°°  and 


q2  -  1 


Properties  (2)  and  (3)  are  sufficient  conditions  for 

2 

the  existence  of  one  resonant  frequency,  for  any  K  ,  provided 
2  2 

that  >1.  Note  that  =1  corresponds  to  the  case 

where  the  primary  circuit  is  non-resonant,  as  was  discussed 
in  Sec.  6.1. 

Properties  (2) ,  (5) ,  (6) ,  and  (7)  are  sufficient 

2 

conditions  (if  Q-.  >  1)  for  the  existence  of  one  resonant 

2  °2 

frequency  in  the  region  0  <  u  <  q  ,  provided  that 


•  ' 


■ 


a93o3qm\**  ow t  E£ri  (  “  '  ' 


& 


evii.^n  ax  **  rloxrlw  needed  <8fr.»  -P*  # 


■  ,b  c-  •  /  Y«  »*  ,*>«.«*•**  S"°  '°  3 

,  ,E  ::  ub  31. w  a*  «  -n:  noesn-  o  .  i  ' 
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Property  (6)  is  a  sufficient  condition  that  no 
resonant  frequency  can  exist  in  the  region  between  the  two 
as^ymptotes . 

Properties  (6)  and  (7)  are  sufficient  conditions 
2 

that,  for  K  =  1 ,  no  resonant  frequency  can  exist  for  which 


.  .  2 

Since  for  every  finite  value  of  u  beyond  the  second 
2  2 

as/fymptote,  F  (u  ,K  )  is  positive  and  finite,  properties  (3) 

and  (6)  are  sufficient  conditions  for  the  existence  of  two 

2 

resonant  frequencies  in  this  region,  providing  Q1  is  suf¬ 
ficiently  large. 

2  2 

A  sketch  of  F (u  ,K  ),  corresponding  to  the  range 


2 

Q. 


2 

<  K  <  1,  is  shown  in  Fig.  6.4.  The  portions  of  the 


curve  represented  by  dashed  lines  indicate  that  the  exact 
nature  of  the  function  in  these  regions  is  not  known.  It  is 
only  known  that  the  function  must  be  continuous  and  of  fixed 
sign,  so  that  the  possibility  of  relative  extrema  has  not  been 
excluded . 


2  1 

2  2  2  . 
The  exact  nature  of  F(u  ,K  )  for  K  <  —  -  — j  /  is 

Qo  Qo 


very  difficult  to  determine.  If  the  function  is  differentiated 


Si 


_  _  r  xi  h  .51091 


-  '  v«  .  n'-'* 
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Figure  6.4 

2  2 

Approximate  Form  of  F(u  ,K  ) 

and  equated  to  zero,  the  resulting  equation  is  considerably 

2 

more  complex  than  the  original  cubic  in  co  (Eq.  6.6).  Thus, 

such  an  undertaking  is  impractical. 

Assuming  as  50  and  Qg  as  /50  ,  F(u  ,K  )  has 

been  evaluated  by  means  of  a  digital  computer,  and  the  re- 

2 

suiting  curves  (corresponding  to  different  values  of  K  )  are 

shown  in  Fig.  6.5.  Only  positive  values  of  the  function  are 

2 

shown,  and  only  positive  values  of  u  are  considered.  The 
curves  are  seen  to  agree  qualitatively  with  the  sketch  in 
Fig .  6.4. 


~  , ,  .  sr  1  '  h  ’■ 
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Note  that  for  Q0  =  50,  =  —  -  -  corresponds  to 

Q,  q7 

2  9  2  2 
K  =  0.0396.  Thus  the  curve  for  K“  =  0.01  should  be  con¬ 
tinuous,  and  all  others  should  have  two  vertical  as/ymptotes. 
This  is  seen  to  be  the  case. 

The  reason  for  choosing  considerably  larger  than 
Qe  is  discussed  in  Sec.  6.7. 

Careful  inspection  of  Fig.  6.5  indicates  that  two 

2 

uniquely  different  situations  can  occur  as  K  is  increased 

2 

from  zero,  depending  on  the  value  chosen  for  .  Fig.  6.6(a) 

2  2  2 

and  (b)  illustrate  the  nature  of  F(u  ,K  )  for  very  small  K  , 
as  can  be  seen  in  Fig.  6.5. 


Figure  6 . 6 

2  2  2 
Approximate  Form  of  F(u  ,K  )  for  Very  Small  K 


oi  abnoq8»Tnoo 
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Note  that  the  relative  maximum  appears  to  occur  for  u2  slightly 
less  than  one.  This  is  approximately  the  same  place  as  the 
as/ymptotes  first  appear. 

Consider  first  Fig.  6.6(a) .  With  reference  to 

•  •  9 

Fig.  6.5,  this  could  correspond  to  -  80.  It  is  apparent 

2  2 

that  for  K  =  0,  u  (for  resonance)  is  considerably  larger 

than  one.  Since  for  K2  =  0,  F(u2,K^)  =  1  +  Q  2u2,  this  value 
2 

of  u  is  determined  by  the  equation 


1  + 


„  2  2 
Q  u 
e 


2 

or  u 


1 


1 

—)  .  (From  6.31) 


2  2  2  2 
If  >>  1,  u  -  Ff  .  But  since  u  is  considerably  larger 

2 

than  one,  this  implies  that  F^  >1. 

2 

Now,  as  K  is  increased,  it  is  apparent  that  the 

original  resonant  frequency  will  be  continuously  increased, 

2 

and  that  for  some  value  of  K  ,  two  additional  resonant  fre- 

2 

quencies  will  appear  at  u  -  1.  Upon  subsequent  increase  of 
2 

K  ,  these  two  frequencies  will  be  separated  by  the  two  as^ymp- 

totes,  as  shown  in  Fig.  6.4.  Thus,  the  loci  of  the  three 

2 

resonant  frequencies,  if  plotted  against  K  ,  will  look  approxi¬ 
mately  as  shown  in  Fig.  6.7(a).  The  three  loci  have  been  de¬ 
signated  as  State  (1) ,  State  (2) ,  and  State  (3) . 


»  KS  -  fl  .  f  <« 


.  .  ■  r 


u 


x  -  3  43  esilqi  i  BiriJ  , -••  '  : 
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Figure  6.7 

Resonant  Frequency  Loci 

Now  consider  the  second  situation,  as  shown  in 

2 

Fig.  6.6(b).  This  could  correspond  to  Q  -  40.  In  this 
2  2 

case,  for  K  =  0,  u  is  seen  to  be  less  than  one.  Again, 

for  =  0  ,  u^  -  F^  (if  q^2  >>  l)  ,  so  that  ¥  ^  must  clearly 

be  less  than  one. 

2 

As  K  is  now  increased,  the  original  resonant  fre¬ 
quency  is  decreased,  and  two  additional  resonances  will  be 

2 

introduced  for  u  somewhat  greater  than  one.  Again,  the  two 
lower  frequencies  will  be  separated  by  the  two  as^ymptotes 


A  _  V  ffv*  .A3£3 
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(Fig.  6.4);  so  that  the  three  loci  can  be  expected  to  look  as 
shown  in  Fig.  6.7(b). 

Since  the  exact  locations  of  the  relative  maximum 

and  minimum,  (Fig.  6.6)  when  they  first  appear,  are  not  known, 

the  exact  border-line  between  the  two  cases  illustrated  in 

Fig.  6.7  cannot  be  determined.  However,  if  F(u  ,K")  is  evalu- 
2 

ated  at  u  =1,  the  result  is 

F(1,K2)  =  Q~2  (  — —  +  K2Q  )2  +  Q2  (6.40) 

2  0  6  6 
u2 

2  2  2 
For  K  =0,  this  lies  on  the  straight  line  1  +  u  ,  which 

of  course  it  must. 

2  2 
Now,  as  K  is  increased,  (6.40)  shows  that  at  u  =  1, 

2  2 

F (u  , K  )  is  continuously  increased.  This  indicates  that  the 

2 

relative  minimum  (Fig.  6.6)  must  lie  to  the  right  of  u  =1. 

2  2 

Furthermore,  a  solution,  for  K  =  0 ,  at  u  =1,  must  be 

2 

initially  pushed  to  the  left  as  K  is  increased,  and  must 
therefore  correspond  to  the  situation  shown  in  Fig.  6.7(b). 

For  such  a  solution, 


or 


1  +  Q. 


Q1  ;  (from  Eq .  6.31) 


1  + 


,  from  which 


Q- 


1 


(6.41) 


■ 
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Thus,  it  can  be  concluded  that  if  Ff* 1 2  <_ 


Ql  -  1 


,  the  resonant 


frequency  loci  must  correspond  to  those  shown  in  Fig.  6.7(b); 

and  for  the  loci  shown  in  Fig.  6.7(a),  F,2  must  be  somewhat 

9  f 

larger  than  — ~ ^ -  . 

Qi  -  1 


It  is  now  noticed  that  State  (1)  and  State  (3)  of 

Fig.  6.7(a)  qualitatively  resemble  the  experimentally  observed 

pattern  shown  in  Fig.  1.5(a).  Also,  the  above  conclusions  con- 

2 

cerning  the  magnitude  of  F^  agree  with  the  experimental  ob¬ 
servation  mentioned  in  Sec.  1.4.  It  will  in  fact  be  shown 
(in  Sec.  6.5)  that  the  above  situation  can  give  rise  to  the 
switching  phenomenon  discussed  in  Sec.  1.4.  However,  before 
this  can  be  done,  the  resonant  resistance  loci,  corresponding 
to  States  (1) ,  (2) ,  and  (3)  must  be  determined. 


6.3  Resistance  at  Resonance 


Assuming  that  a  numerical  (or  graphical)  solution  to 
Eq.  6.30  has  been  found,  the  impedance  at  the  resulting  reso¬ 
nant  frequencies  must  now  be  determined. 

R  '  shall  be  defined  as  the  resistance  of  any 
o 

2 

resonant  state,  as  determined  by  Eq.  6.30.  Thus,  for  K  —  0, 

R  '  =  R  .  Also,  G  '  shall  be  defined  as  1/R  '  . 
o  o  o  o 

Referring  now  to  Eq.  6.25,  at  resonance 


YT(  jw) 


(R,  +  R  ) 

1  0) 


2  2  2 

(R,  +  R  ;  +  w  L 

1  00  00 


I 


(6.42) 


■ 
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Also ,  from  (6.26)  , 


00 


2  2  2' 
( R,  +  R  )  +  id  L 

1  o>  co 


which  can 


be  rearranged  as 


2  2 

co  L 

oo 


-  -  <Ri  +  V2 


(6.43) 


Substituting  (6.-43)  into  (6.42), 


G  ' 

o 


(Rn  +  R  ) 

1  00 


(R1  +  RJ2  +  57 


(Rl  +  V 


(R,  +  R  )C. 

1  00  1 


or 


00 


R  ' 

o 


00 


(R1  +  RU)C 


(6.44) 


Recalling  the  definition  of  and  R^  (Eq's.  6.17 
and  6.18),  (6.44)  can  be  written  as 


R  ' 

o 


LnL 
1  n 


C1R1(1  +  Rn> 


But 


R1CL 


R  ,  so  that 
o 


R  '  =  R 

o  o 


n 


1  +  R 


or 


n 


R  ' 

o 


n 


R 


o 


1  +  R 

n 


(6.45) 


■ 
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Substituting  (6.19)  and  (6.24)  into  (6.45), 


R  ' 
o 


(1  -  u  ) 


1  -  u* 2 (1  -  K2)|  +  u 


■] 


Qo2 


(6.46) 


R 


o 


(1  -  u2,2  + 


u 


4  2  Qe 
+  u  K  — 


Q. 


It  must  be  kept  in  mind  that  Eq.  6.46  is  valid  only 
2 

for  values  of  u  as  determined  by  Eq.  6.30. 


6.4  Discussion  of  the  Function  R  ' /R 

o  o 


It  should  first  be  pointed  out  that  only  the  same  two 

2  2 

parameters,  Q?  and  Q  ,  appear  in  this  function  as  in  F(u  ,K  ). 
This  is  very  important  with  respect  to  optimization,  as  it 
implies  that  both  functions  can  be  evaluated  for  a  specified 
Qe  and  Q^,  Furthermore,  if  (and  F^)  are  varied,  the  cor¬ 
responding  variations,  both  in  resonant  frecruencies  and  in 


R  ' 
o 


resonant  resistances,  can  be  studied  simultaneously. 

It  can  easily  be  seen  from  (6.46)  that  lim 

K2+0  o 

This  of  course  corresponds  to  a  single  tank  circuit.  Also 


=  1 


lim 

u^+0  R 


R  1 

o 


=  1  ,  and 


(6.47) 


o 


it 


im 


R  ' 
o 


i  -  r 


u  ->°°  R 


o  Q 

1  +  K  Ql 


(6.48) 


■  IB- 
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Equation  6.48  represents  a  horizontal  as^ymptote,  which  must 

2 

lie  between  zero  and  one,  for  0  <  K  <1. 

Since  the  denominator  of  Ro'/R0  must  always  be  posi¬ 
tive,  and  the  numerator  is  the  second  factor  in  the  denominator 

2  2  2  2 
of  F(u  , K  ) ,  Rq ' /R  must  have  the  same  sign  as  F(u  , K  ) ,  and 

2  2 

must  be  zero  where  F(u  ,K  )  has  vertical  as^ymptotes.  Thus, 
at  any  solution  to  Eq.  6.30,  Rq ' /Rq  must  be  greater  than  zero. 

If  Rq ' /R  is  equated  to  one,  the  resulting  equation 


is 


(1-U2) 


2  2 
l-uz(l-Kz) 


u 


+ 


2  2 

=  (1-u  )  + 


u 


.  4  2  Qe 

+  u  K  - 


Q. 


Upon  simplification,  this  can  be  expressed  as 


u 


1  -  u2(l  +  — ) 

Q 


2  -1 


=  0;  for  K  ^  0 


Thus , 


u  =  0 


or 


2 

u  = 


Q, 


(6.49) 


1  + 


Since  there  are  only  two  points  for  which  Ro'/jR0  =  1, 


and  since 


i1 


R  ' 

o 


1  -  K‘ 


.im 

U“->oo  R 


2  Qe 
1  +  K  — 


,  it  follows  that 


9  ■  ;  nojJi-iipa 


■ 

a- 


if 


0  <  u  < 


1  + 


5 

Q. 


-  94  - 


R  ' 

,  then  —  ■  >  1 ;  and 
R 


(6,50) 


if 


2 

u  > 


1  + 


Q. 


R  ' 

,  then  -2_  <  io 


R 


o 


Referring  now  to  Fig.  6.4,  (6.50)  implies  that  for  the  two  pos 

sible  solutions  to  the  right  of  the  second  asymptote,  Ro'/R0 
must  be  less  than  one,  and  for  the  solution  to  the  left  of  the 
first  asymptote,  Rq’/Ro  can  be  either  greater  or  less  than  one 

It  is  also  of  interest  to  evaluate  the  function 
2 

R0'/R0  at  u  =1.  The  result  is  seen  to  be 


R  ' 

o 

R 


2  n 

u  =1 


1  +  k  QeQ2 


(6.51) 


Figure  6.8  illustrates  the  general  nature  of  RoVR0 


2  2 

(as  concluded  from  the  previous  discussion)  for  K  >  — 

Q. 


2  ° 


2  2 

The  relationship  to  the  function  F (u  ,  K  )  is  also  shown.  The 
dashed  lines  again  indicate  regions  of  relative  uncertainty 
as  to  the  exact  nature  of  the  functions. 

It  is  apparent  from  Fig.  6.8,  that  Rq Vr0  must  have 


2 

an  absolute  maximum  in  the  region  0  <  u  < 


1  +  V°2 


and 


£  I 
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Figure  6.8 

Approximate  Form  of  Rq ' /R  and 

2  2 

its  Relationship  to  F(u  ,K  ) 


an  absolute  minimum  in  the  region  where  the  function  is 

negative  (for  K2  >  — - .  The  ordinate  values  of  these 

Q2  Q2  2 

extrema  can  quite  easily  be  found  for  K  =1.  This  will  be 
of  further  assistance  from  the  point  of  view  of  optimization. 


Letting  K  = 


(l-u2)  +  U 


R  ' 
o 


R 


o 


Q. 


K2=l 


2  2  u‘ 
(1-uV  +  - 


u 


4Q, 


+ 


lr 
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which  can  be  rearranged  as 


R  ' 
o 

R 


1  -  U4(l - t) 


(6.52) 


K2  =  l 


4  ^e  2  ^ 

u  (1  +  — )  -  u  (2 - tt)  +  1 


^2  v2 

2 

If  x  is  substituted  for  u  in  Eq.  6.52,  the  resulting  equation 
is  of  the  form 


1  -  ax 

y(x)  =  — « -  (6.53) 

bx^  -  cx  +  1 

where  a,  b,  and  c  are  implicitly  defined  by  (6.53)  and  (6.52). 

If  Eq.  6.53  is  equated  to  an  arbitrary  constant,  X,  the  re- 

2 

suiting  equation  will  be  (1  -  ax)  =  A (bx  -  cx  +  1) ;  which  can 
be  expressed  as 

x2Xb  -  x(Ac  -  a)  +  (A  -  1)  =  0  (6.54) 

If  A  is  now  taken  to  represent  an  extreme  value  of  R  ' /R  , 
there  can  be  only  one  solution  to  the  equation 


x  Ab  -  x ( Ac  -  a) 
(Ac  -  a) 2  = 


+  (A 
4  Ab  ( A 


1)  = 
1)  . 


0;  from  which 

This  can  be  rearranged 


as 

A2 (c2  -  4b)  +  A(4b  -  2ca)  +  a2  =  0  .  (6.55) 

The  two  roots  to  (6.55)  are 

-(4b  -  2ca)  ±  /(4b  -  2ca) 2  -  4a2 (c2  -  4b) 

A  - - ^ 

2(c  -  4b) 
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which  simplify  to 

(ca  -  2b)  +  2/b?  -  abc  +  a2b 


A  = 


-  4  ab 


( 6  o  56 ) 


If  the  values  of  a,  b,  and  c  are  now  substituted 
into  Eq.  6  0  56  an  exact  expression  for  A  can  be  obtained.  How¬ 


ever,  if  Q_  >>  1,  then  a  -  1,  and  c 


No  such  approxi¬ 


mation  can  be  made  for  b,  so  that  b  =  1  +  Qe/Q2 °  Substituting 
these  values  into  ( 6  0  5 6 )  and  simplifying, 


A 


-  (1  +  /I  +  Q,/Q J 
2  *  e 


( 6  o  57 ) 


Equation  6.57  indicates  that  the  maximum  possible 
resonant  resistance  is  determined  primarily  by  the  ratio 


W 

The  function  R  ' /R  has  been  evaluated  for  Q„  =  50 

o  o  2 

and  Qe  =  /50*,  and  is  shown  in  Fig.  6.9.  Substituting  the 

above  values  into  (6.57),  A  -  1.92  and  -0o92o  These  are 

2 

seen  to  agree  with  the  peaks  shown  in  Fig.  6.9  (for  K  =1) 0 

The  resonant  resistance  loci,  corresponding  to  the 

three  resonant  states  shown  in  Fig.  6.7(a),  can  now  be 

qualitatively  determined.  Consider  first  the  locus  designated 

2 

as  State  (1) .  Since  for  this  state,  u  is  always  greater  than 

2 

one,  it  is  apparent  from  Fig.  6.9,  that  for  K  >0,  Ro8/R0  < 

2  2  2  i  _  j^2 

Also,  as  K  ->  1,  u  -*  and  as  u  -*  °°,  R  */R  ->  - — - — . 

1+K2  Qe/Q2 

2 

This  approaches  zero  as  K  approaches  one.  Thus,  the  locus 


»  ;  +  r:t8 dus  -won  c o  fens  ,d  *B  *0  v  91 'l*  '  1 


.pnx'flii^raiB  bne  Oa.d).  o3ni  bsi/Xbv  assrl^ 


bns  se.X  «  A  ,(«.d)  oinx  bsuXbv  avode 
„(X»:X  ioi)  e  a  .pi's  ni  nworia  8>[B9q  £*<ii  ri:*iw  ©sips  oi  r 
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Figure  6 . 9 
Graph  of  R^/Rq 
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2 

°f  R0'/R0'  plotted  against  K  ,  must  look  somewhat  as  shown  in 
Fig.  6.10,  State  (1)  . 


Figure  6.10 

Approximate  Resonant  Resistant  Loci 

Now  consider  States  (2)  and  (3),  of  Fig.  6.7(a). 

These  two  states  are  introduced  simultaneously  at  some  value 
2  2  2 

of  u  .  If  (and  F^  )  are  such  that,  at  this  point, 

,  it  can  be  seen  from  Fig.  6.8  that  R  ' /R 


2 

u  > 


1  +  V°2 


must  initially  be  less  than  one.  Also,  for  State  (2),  as  K  +1 , 
2 

u  ,  so  that  Rq  1  /Rq  must  approach  zero,  as  for  State  (1)  . 

It  can  be  seen  from  Fig.  6.8,  that  if 


1  2 

— 2  <  K  <  1/  State  (2)  will  be  closer  to  the  second 


^2  ^2 

as^ymptote  than  State  (1) ,  and  thus  will  have  a  lower  resonant 

qT  '  o7 


2  2 

resistance.  This  will  also  be  true  for  K  <  — 


if 


State  (2)  is  near  the  minimum  of  Ro'/R0  (Fig.  6.9). 


l'  1  ' 


. 


■ 
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Thus,  the  locus  of  Rq '/Rq  for  State  (2)  must  be  approximately 


as  shown  in  Fig.  6.10. 


2  2 

For  State  (3),  as  K*  increases,  u 


2  2 

if  for  some  value  of  K  ,  u  < 


R 


1  +  Qe/Q2  '  o 

greater  than  one.  The  final  value  of  R  ‘/R 

.  o  '  o 

be  determined  by  the  point  on  the  curve 


decreases.  Hence, 

/R  must  be 
'  o 

2 

as  K  +1,  will 


R  ' 

o 


2 

corresponding  to  the  final  value  of  u  for  State  (3) .  From 

Eq.  6.57,  this  must  be  less  than,  or  approximately  equal  to 

%(1  +  /l  +  ) •  The  approximate  shape  of  the  locus 

R  ' /R  for  State  (3)  is  shown  in  Fig.  6.10. 

By  careful  inspection  of  Fig's.  6.5  and  6.9,  it 

can  be  seen  that  the  general  forms  of  loci  indicated  in 

2 

Fig's.  6.7(a)  and  6.10  can  be  realized  if  is  chosen  as 

2  2  2  2 

70.  Since  F^  =  /Qe  ,  and  Qg  =50,  this  implies  that 

Ff  2  =  70/50  =  1.4  . 

Using  Q2  =  50,  Qg  =  /50",  and  Q1  =  /70,  a  digital 
computer  has  been  used  to  numerically  solve  Eq.  6.30  (for 
different  values  of  K  ) ,  and  subsequently  to  evaluate  R  ' /R  . 
The  resulting  resonant  frequency  and  resonant  resistance  loci 
are  shown  in  Fig's.  6.11  and  6.12.  Note  that  RQ ' /RQ  is  much 
less  for  State  (2)  than  for  State  (3) . 
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Resonant  Frequency  Loci 
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Resonant  Resistance  Loci 
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6 • 5  Explanation  of  Switching  Phenomenon 


The  previous  theory  concerning  the  resonant  states 
for  the  mutually  tuned  circuits  can  now  be  applied  in  con¬ 
junction  with  the  describing  function  method  as  developed  in 
Sec's.  4  and  5.  An  explanation  of  the  switching  phenomenon, 
discussed  in  Sec.  1.4  will  thus  result. 

For  the  oscillator  circuit  shown  in  Fig.  1.4,  the 
basic  equations  for  sinusoidal  oscillation  will  now  be 


G 

eq 


<A,V  ) 


<A'V 


1 

- -  ,  and 

R  ’ 

o 


+  R 


1 


(6.38) 


Thus,  a  solution  in  A  and  will  be  determined  in  exactly 
the  same  manner  as  described  in  Sec.  5,  except  that  the  con¬ 
tour  G  =  -  1/R  '  must  be  chosen, 
eq  '  o 

The  frequency  of  oscillation  will  be  given  by 
2  2  2 

co  =  u  /L~C~  ;  where  u  is  a  solution  to  Eq.  6.30. 

o  '  2  2  ^ 

Fiq.  6.13  shows  four  contours  on  a  Constant-G 
r  eq 

Chart,  and  the  curve  Vq(A)  for  some  load-line.  This  could 
correspond  to  (R-j^  +  Rg)  =  1250,  and  Vg  =  237.5  mv,  as  shown 
in  Fig's.  5.7  and  5.8.  Also,  Fig.  6.14  illustrates  the 
resonant  frequency  and  resonant  resistance  loci  of  Fig's.  6.11 


and  6.12. 


. 


I  b  i  S  3  wen  >  i  4  u  '  ■' 


x 


: 


104 


eq 


Resonant  Frequency  and  Resonant  Resistance  Loci 
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With  reference  to  Fig.  6.13,  assume  that  -  1/R  =  r3.  Thus, 

2 

for  K  =0,  oscillation  will  occur  with  amplitude  A^ .  Also 

in  Fig.  6.14(b)  assume  that  at  and  P^,  -  1/R  '  =  r^;  at 

P2,  -  1/R0'  =  r4;  and  at  P4,  -  1/RQ '  =  r2- 

2  . 

Now,  as  K  is  initially  increased,  Rq  '  for  State  (1) 

will  decrease  until  P1  is  reached  (Fig.  6.14(B)).  This  will 

correspond  to  an  increase  in  frequency,  as  indicated  by  Fig. 

6.14(a).  Also,  the  solution  to  Eq's.  6.58  will  move  along 

the  curve  Vq (A) ,  until  the  countour  is  reached;  which  will 

correspond  to  oscillation  of  amplitude  A^ .  It  is  apparent 

2 

that  upon  subsequent  increase  of  K  ,  no  solution  to  Eq's. 

6.58  can  exist  for  the  value  of  R  '  corresponding  to  State  (1) . 

2 

However,  for  this  value  of  K  ,  -1/R  ' ,  for  State  (3) ,  is 
equal  to  r^,  so  that  oscillation  may  take  place  with  ampli¬ 
tude  A^ ;  and  frequency  as  determined  by  P2,  Fig.  6.14(a). 

2  2 

Thus,  at  K  =  K2  ,  the  oscillator  can  be  expected  to  switch 

2 

modes.  Further  increase  in  K  will  now  result  in  a  decrease 

in  frequency,  and  an  increase  in  amplitude. 

If  is  now  decreased,  R  '  for  State  (3)  will 

decrease  until  P3  is  reached.  This  will  correspond  to  an 

increase  in  frequency,  as  indicated  by  Fig.  6.14(a);  and  a 

2  . 

decrease  in  amplitude,  until  A^  is  again  reached.  If  K  is 
decreased  beyond  this  point,  no  solution  to  Eq's.  6.58  for 
State  (3)  can  exist.  However,  -  1/RQ '  for  State  (1)  is  now 
equal  to  r2,  so  that  oscillation  of  amplitude  A2  may  occur. 


. 
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Thus,  the  oscillator  can  be  expected  to  switch  back  to  its 

.  .  2  2 
original  mode  at  K  =  K  * . 

It  should  be  pointed  out  that  for  State  (2),  R  ' 
is  always  less  than  the  minimum  value  for  sustained  oscilla¬ 
tion.  (i.e.  the  value  at  P^  and  P  ,)  .  Thus,  oscillation  in 
State  (2)  can  never  occur. 

Example  6.1 

Fig's.  5.7,  5.8,  and  6.12  can  now  be  used  to  design 
a  system  which  should  exhibit  the  switching  pheno¬ 
menon.  Referring  to  Fig.  5.7,  +  R^  can  be 

chosen  as  125ft,  and  as  237.5  mv,  so  that  all 

a 

solutions  in  A  and  will  be  on  the  fourth  lowest 
curve,  Vq(A).  (Fig.  5.8). 

3 

If  Rq  is  chosen  as  1/1.8(10  )ft,  so  that  -1/Rq  = 

2 

-1.8  mmho,  then  for  K  =0,  the  circuit  can  be 
expected  to  oscillate  with  A  =  177.5  mv  and  = 

185  mv. 

It  can  be  seen  from  Fig.  5.8  that  the  minimum 

value  of  R  1  for  sustained  oscillation  corresponds 
o 

3 

to  G  =  -4.5  mmho,  for  which  R  '  =  1/4.5(10  )ft. 
eq  o 

Hence,  R  '/R  =  1.8/4. 5  =  0.4.  Referring  to  Fig. 

o  o 

2 

6.12,  for  State  (1),  this  occurs  for  K  -  0.37. 

From  Fig.  5.8,  the  minimum  amplitude  is  seen  to  be 
approximately  82.5  mv . 

2 

Now  from  Fig.  6.12,  State  (3) ,  for  K  = 


0.37, 


■ 
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R0'/R0  =  1.44,  so  that  -1/RQ'  =  -1.8/1.44  mmho  - 

-1.25  mmho.  Thus,  from  Fig.  5.8,  oscillation  in 

State  (3)  may  occur,  with  A  -  195  mv  and 

VQ  =  192.5  mv. 

2 

If  K  is  now  decreased,  oscillation  in  State  (3) 

will  continue  until  R  '/R  for  State  (3)  reaches 

o  o 

0.4,  so  that  the  minimum  amplitude  (82.5  mv)  is 

again  reached.  From  Fig.  6.12,  this  corresponds 
2 

to  K  -  0.011.  At  this  point,  Ro'/R0  for  State 
(1)  is  approximately  0.95  so  that 
1  -1.8 

-  -  -  -  mmho  =  -1.9  mmho. 

R  '  0.95 

o 

Hence,  from  Fig.  5.8,  State  (1)  oscillation  may 

occur  with  A  -  172  mv  and  -  184  mv. 

The  normalized  frequencies  at  the  switching  points 

2 

can  now  be  found  from  Fig.  6.11.  For  K  =  0.37, 

u2  for  State  (1)  is  about  3.0,  and  u2  for  State (3) 

is  about  0.715. 

2  2 

For  K  =  0.011,  u  for  State  (1)  is  about  1.43, 

2 

and  u  for  State  (3)  is  about  0.98.  The  actual  fre- 

2 

quencies  corresponding  to  the  above  values  of  u 
can  be  found  from  Eq.  6.62. 

It  is  of  interest  to  consider  what  would  happen  if 
the  oscillator  discussed  in  Example  6.1  were  stopped  at  a 
time  when  K  were  somewhere  between  the  two  switching  points. 


. 


'  v- 
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It  can  be  seen  from  Fig.  5.8  that  for  any  contour  for  which 

Geq  <  ~3.0  mmho ,  there  will  be  either  two  solutions  (which 

may  coincide) or  no  solution.  It  has  been  shown  in  Sec.  5.7 

that  when  there  are  two  solutions,  the  one  occurring  for  the 

smallest  amplitude  is  unstable,  and  the  other  one  is  stable. 

Now  for  -  1/Rq '  =  -3.0  mmho,  R  '  =  333ft.  Therefore,  whenever 

V  <  333ft,  for  either  State  (1)  or  State  (3) ,  that  state  will 

be  characterized  by  an  unstable  and  a  stable  solution  (or  no 

solution  at  all).  For  R  '  =  333ft,  R  WR  =  1. 8/3.0  =  0.6. 

o  o  o 

Referring  to  Fig.  6.11,  for  State  (1)  this  corresponds  to 

2  2 
K  =  0.17,  and  for  State  (3) ,  it  corresponds  to  K  -  0.02. 

Thus,  if  0.02  <  <  0.17,  R  ' /R  for  both  State  (1)  and  State 

(3)  will  be  greater  than  0.6,  and  only  stable  solutions  will 

exist.  This  implies  that  if  the  oscillator  were  stopped  and 

permitted  to  start  again,  the  mode  of  oscillation  could  not 

be  predicted. 

2 

If  K  >  0.17,  there  will  be  only  one  (stable)  solu¬ 
tion  for  State  (3) ,  and  an  unstable  and  a  stable  solution  for 
State  (1)  (or  no  solution).  Hence,  in  this  case,  if  the 

oscillator  were  allowed  to  start,  oscillations  would  occur  in 

2 

State  (3).  Similarly,  if  K  <  0.02,  oscillations  would 
spontaneously  occur  in  State  (1) . 

It  is  conceivable,  that  under  certain  conditions, 
oscillation  in  both  modes  may  occur  simultaneously.  If  the 
two  frequencies  were  fairly  close  together,  a  beat-frequency 


. 


-- 
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phenomenon  would  thus  be  exhibited.  Analysis  of  the  system 
under  such  circumstances  would  require  the  simultaneous  solu¬ 
tion  of  three  equations,  in  V^,  A^  ,  and  A^;  where  A^  and  A 2 
represent  the  fundamental  amplitudes  of  the  two  voltage  com¬ 
ponents.  Such  an  analysis  will  not  be  undertaken  in  this 
thesis,  as  it  would  be  much  too  complex  to  be  feasible. 

6 . 6  Limitations  of  the  Method  of  Analysis 

Because  the  describing -function  method  is  valid  only 
for  steady-state  sinusoidal  oscillations,  no  information  has 
been  obtained  regarding  transient  response.  Also,  nothing 
has  been  learned  about  the  exact  nature  of  the  switch  from 
one  mode  to  another.  However,  it  seems  logical  to  assume 
that  transients  of  some  form  would  occur  during  the  switch¬ 
ing  process. 

It  has  also  been  observed  that  some  very  non- 
sinusoidal  modes  of  oscillation  can  occur,  particularly  when 
both  tuned  circuits  are  used.  The  describing-function  method 
is  of  course  inadequate  for  analysis  of  such  modes. 

In  order  for  the  assumption  of  sinusoidal  oscilla¬ 
tion  to  be  nearly  met,  it  has  been  shown  (in  Sec.  5.9)  that 
the  impedance  to  the  harmonics  of  the  frequency  of  oscilla¬ 
tion  should  be  as  small  as  possible  compared  to  rd  for  the 
tunnel  diode.  For  one  tuned  circuit  alone,  this  is  not  dif¬ 
ficult  to  achieve.  However,  when  two  coupled  circuits  are 


. 
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used,  the  frequency  response  of  the  impedance  function  is 
known  to  vary  drastically  as  coupling  is  varied.  Thus,  under 
certain  circumstances,  it  is  to  be  expected  that  varying  the 
degree  of  coupling  will  vary  the  distortion,  and  thus,  the 

frequency  departure.  As  an  example,  if  the  circuit  were 

.  .  2 
oscillating  in  State  (3),  (Fig.  6.12),  and  K  were  somewhat 

larger  than  the  value  required  for  the  circuit  to  switch  to 
State  (1) ,  the  resonant  resistance  for  State  (1)  would  be 
larger  than  that  for  State  (3) .  If  the  resonant  frequency 
of  State  (1)  were  approximately  twice  that  for  State  (3) , 
the  impedance  to  the  second  harmonic  would  be  considerably 
larger  than  the  impedance  to  the  fundamental.  Thus,  consider¬ 
able  distortion  would  probably  result. 

6 . 7  Summary 

The  switching  phenomenon,  as  predicted  by  the  de¬ 
scribing-function  method,  is  seen  to  depend  on  a  combination 
of  circumstances  involving  both  the  linear  and  the  non-linear 
portions  of  the  system.  Regarding  the  non-linear  portion, 
such  a  switch  can  be  predicted  only  when  conditions  are  such 
that  a  stable  and  an  unstable  solution  in  A  and  VQ  may  occur, 
so  that  there  will  be  a  minimum  amplitude  for  sustained 
oscillation.  Thus,  the  curvature  of  the  contours  G0q  =  r, 
and  of  the  loci  V  (A)  are  of  primary  importance. 


' 


* 
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In  general, two  solutions  may  occur  in  any  bias  range 

for  which  the  contours  G  =  r  are  concave  upward  at  A  =  0 

eg 

(if  vq  >  vj)-  ln  Sec.  5.4,  this  range  has  been  identified  as 
the  range  for  which  g^(v^)  <  0.  Also,  by  Theorem  5.1, 

^3  (Vq)  <  0  in  any  range  of  for  which  dp/dS  >_  0. 

If  the  d.c.  load-line  intersects  the  d.c.  character¬ 
istic  slightly  to  the  right  of  the  inflection  point,  the 
curvatures  of  the  contours  and  of  the  curves  V^(A)  are  both 
conducive  to  the  appearance  of  two  solutions. 

For  the  linear  portion  of  the  system,  the  first 
requirement  for  the  switching  phenomenon  to  occur  is  that 
the  resonant  frequency  loci  for  the  mutually  coupled  circuits 
must  resemble  those  shown  in  Fig.  6.7(a) .  It  has  been  shown 

in  Sec.  6.2  (Eq.  6.41)  that  a  necessary  condition  for  this 

2  2  2 
is  that  Ff  be  somewhat  larger  than  Q1  /(Q^  -1) . 

2 

A  second  requirement  is  that  as  K  is  increased 

(or  decreased)  the  second  mode  must  be  available  before  the 

first  mode  becomes  unavailable.  If  this  were  not  so,  upon 

2 

increasing  (or  decreasing)  K  ,  at  some  point  oscillation  would 
cease  entirely,  and  would  re-occur,  in  the  second  mode  when 
some  other  point  were  reached.  Referring  to  Fig.  5.12,  this 
requirement  means  that  Rq '/Rq  for  State  (1)  should  initially 
decrease  as  slowly  as  possible;  and  Rq'/R0  for  State  (3)  should 
increase  as  rapidly  as  possible.  It  has  been  shown  in  Sec. 


6.4  (Eq.  6.57)  that  for  State  (3),R0'/R0  £  Hd  +  /l  +  Q2/Qe  )  * 


. 
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2 

Also,  from  Fig.  6.8,  it  is  seen  that  for  State  (3),  u  must 
be  less  than  ■= — -<-  ^ — jrr—  ,  in  order  that  R  '  /R  be  greater 

1  +  Qe/Q2  o  '  o  ^ 

than  one.  Thus,  R^'/R^  will  increase  most  rapidly,  and  have 

2 

the  largest  value  for  K  =1,  when  Q9/Q  is  made  as  large  as 

M  6 

possible.  Further  optimization  in  this  respect  can  be  ob- 

2  2 

tamed  by  choosing  Q1  (and  F^  )  so  that  the  solutions  to  Eq. 
6.30,  for  State  (3),  be  approximately  on  the  peaks  shown  in 
Fig.  6.9.  This  is  the  basis  on  which  Q Q^,  and  F^  were 
chosen  for  Fig's.  6.5,  6.9,  6.11,  and  6.12. 
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VII.  EXPERIMENTAL  VERIFICATION 

An  analogue  computer  can  easily  be  used  to  simulate 
a  tunnel  diode  oscillator,  and  thus  provides  a  convenient 
means  of  testing  the  previous  theory.  Since  the  frequency 
of  oscillation  is  irrelevant  to  the  nature  of  the  solution, 
the  computer  simulation  can  be  designed  to  oscillate  at  any 
desirable  frequency.  Thus,  the  time-domain  solution  can  be 
effectively  slowed  down  so  that  transient  response  can  be 
observed.  Furthermore,  the  actual  switching  process  can 
be  studied. 

A  further  advantage  to  the  analogue  computer  simu- 

2 

lation  is  that  the  coupling  coefficient,  K  ,  can  accurately 
be  adjusted  to  any  value  by  adjusting  a  single  pot.  This  is 
shown  in  Sec.  A2 . 4 . 

The  analogue  computer  simulation  can  of  course  be 
used  to  study  non-sinusoidal  modes  as  well,  if  this  be  de¬ 
sired  . 

7 . 1  Basic  Computer  Set-Up 

Consider  the  block-diagram  of  the  hypothetical 
control  system  shown  in  Fig.  7.1.  Two  basic  equations  can 
be  written  for  this  system; 

x2  =  f  (x^  f  and  (7.1) 

=  R  ( S )  -  x 2  (S)  {E±  (S)  +  H2  (S)  } 


x1  (S) 


(7.2) 


•  • Joqyri  io  ti£-  c  tb-jtoold  mU 


' 

iMttW.V*  « «t*i»  9d 
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Figure  7.1 

Block  Diagram  of  Computer  Simulation 

Equations  (7.1)  and  (7.2)  are  identical  to  Eq ' s .  (4.11)  and 

(4.10) ,  if; 

R  ( S )  =  VB(S);  x^  (S)  =  VD(S);  x2(S)  =  ID  (S)  ; 

H1(S)  =  ZT(S);  and  H2(S)  =  ZB(S). 

Thus,  accurate  simulation  of  the  circuit  shown  in  Fig.  4.3 
requires  only  a  non-linear  function  generator  which  will 
produce  the  d.c.  characteristic  for  the  tunnel  diode;  a 
transfer  function  mathematically  equal  to  the  impedance 
function  of  the  tank  circuit;  and  a  transfer  function  mathe¬ 
matically  equal  to  the  impedance  function  of  the  biasing  net¬ 
work.  When  mutual  coupling  is  included,  H^(S)  must  of  course 
be  equal  to  ZT(S),  as  shown  in  Fig.  6.1. 

Detailed  computer  schematics  for  the  above  system 
components  are  derived  in  Sec .  A2 . 


. 
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All  experimental  results  included  in  this  section 


have  been  obtained  using  a  PACE  231-R  analogue  computer  and 
an  1100-E  X-Y  recorder. 

7 . 2  Oscillator  Without  Mutual  Coupling 

The  computer  schematic  shown  in  Fig.  A2 . 8  can  be 
used  to  investigate  any  aspect  of  the  describing-function 
theory,  as  well  as  the  various  non-sinusoidal  modes  of  oscil¬ 
lation  that  may  occur. 


The  circuit  discussed  in  Example  5.1  shall  first 


be  simulated.  For  this  circuit,  R  =  250ft,  and  (R,  +  R_.)  = 

o  -L  a 

125ft.  The  other  circuit  parameters  can  be  chosen  as  desired. 


Choosing 


10 


,  and 


(7.3) 


1 


(7.4) 


2tt/1  -  1/100  ;  or  f  =  .  995  Hz. 


From  Eq ' s .  3.4  and  3 . 7 


R 


o 


or 


/Ll/Cl  =  25ft. 

From  (7.4)  and  (7.5),  C±  =  1/50tt  ,  and  L1  =  100/8tt  . 


(7.5) 


. 


.  A  4 


- 
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Also,  from  Eq.  3.7,  Rl  =  £T“  ^l^l'  from  which  =  2.5ft. 

Since  (R^  +  Rg)  =  125ft,  Rfi  =  122. 5ft.  The  impedance  of  Cg  can 
be  neglected  if 


. 01  R  ,  or 

o 


100  100  1 

^  a)  R  2 tt  x  250  5tt 

o  o 

Note  that  /L^/C^  -  7  x  10  ^  x  25  =  .175,  so  that 
near  sinusoidal  oscillation  can  be  expected  to  occur.  (see 
Sec ' s .  5.9  and  7.3.) 

If  the  current  scale  is  chosen  as  1  ma  =  1  V, 

3 

f (VD)  will  be  multiplied  by  10  ft,  so  that  all  linear  impe- 

3 

dances  must  be  divided  by  10  ft.  Hence,  for  the  computer 
schematic , 


R 


1 


L 


1 


C 


1 


.025 

1 


8  0  7T 

20 

7T 


c 


B 


.1225, 

200 

TT 


and 


Using  the  above  values,  the  tracings  shown  in 
Fig's.  7.2  and  7.3  have  been  obtained.  Fig.  7.2  shows  the 


wJt/qmoo  ®ri 3  ,aonaH  ,8  01  bzblvlb  ad  *aum  eaonab 

'  i  '  ,oiis  !oa 
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Figure  7.2 

Portions  of  Diode  Characteristic  Covered 
During  Steady-State  Oscillation 


■ 
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portion  of  the  diode  characteristic  covered  during  steady- 
state  oscillation,  corresponding  to  the  values  of  V_  as  in- 
dicated  on  the  curves.  A  and  can  be  obtained  directly 
from  these  curves,  and  are  tabulated  in  Table  7.1 


Table  7.1 

Experimental  Results 


B 

A  (mv) 

V  (rav) 

199 

65 

100 

200 

74 

106 

205 

89 

119 

210 

98 

127 

215 

102 

137 

220 

107 

142 

225 

111 

149 

230 

114 

155 

235 

117 

162 

240 

117 

168 

245 

113 

176 

247 

109 

180 

Oscillations  were  seen  to  die  out  when  V_.  was  in- 

B 

creased  to  248  mv  or  decreased  to  198  mv,  so  that  the  minimum 
amplitudes  for  sustained  oscillation  are  approximately  109  mv 
and  65  mv.  The  above  results  agree  with  the  predictions  of 
the  describing-function  method.  (See  Example  5.1  and  5.2) 


■ 


i 


mmluim  «mW  Jwtt  ob  ,  vm  8CX  at  **8B«soeb  10  Vnl  8"'s  °*  b  E  19 '  - 
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Typical  Solutions 
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Figure  7.3  shows  three  typical  time-domain  solu¬ 
tions  for  VD .  For  Fig.  7.3(a),  Vg  has  been  set  at  225  mv, 
so  that  oscillations  built  up  until  the  stable  solution  was 
reached.  can  be  seen  to  decrease  slightly  as  A  increases. 

In  Fig.  7.3(b),  the  system  was  initially  oscillat¬ 
ing,  and  Vg  was  then  increased  to  250  mv,  so  that  the  solu¬ 
tion  disappeared.  Thus,  oscillations  died  out.  As  the  ampli¬ 
tude  decreases  is  seen  to  increase,  as  predicted  by  the 
describing-function  theory.  (See  Example  5.2) 

In  Fig.  7.3(c)  Vn  has  been  decreased  to  190  mv, 

a 

so  that  again  the  solution  disappeared.  In  this  case,  as 
amplitude  decreases  also  decreases. 

The  steady-state  portion  of  Fig.  7.3(a)  is  very 
nearly  sinusoidal  and  the  frequency  of  oscillation  is  approxi¬ 
mately  .995  Hz  as  was  predicted.  The  final  amplitude  is  about 
112  mv. 

7 . 3  Distortion  With  an  Infinite-Q  Circuit 

By  letting  R^  =  0  for  the  computer  simulation, 

the  dependence  of  distortion  and  freauency  departure  on  the 

quantity  G^/L^/C^  can  be  studied  (see  Sec.  5.9).  For  =  0 , 

both  Q,  and  R  will  be  infinite,  so  that  for  sinusoidal 
1  o 

oscillation,  all  solutions  should  be  on  the  contour  Geg(A,VQ) 

=  0;  and  frequency  should  be  such  that  coo  =  1//L^C^. 


c  V. 


>fo  o  if  r  Qv  898  5010!  /.  sbul  1 J  q^lfi 


■ 
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Also,  in  this  case,  the  biasing  network  can  be  neglected,  so 
that  =  VD. 

In  varying  G^/L^/C^,  it  is  desirable  that  u>  remain 
constant,  so  that  the  frequency  departure  can  easily  be  ob¬ 
served.  This  can  be  conveniently  achieved  by  multiplying 
f (Vp)  by  different  values,  which  corresponds  to  varying  the 
gain  in  the  forward  loop  of  the  computer  simulation. 

Using  L,  =  1/80tt,  C,  =  20/tt,  and  =  V  =  250  mv, 
the  tracings  shown  in  Fig.  7.4  have  been  obtained.  The  cur¬ 
rent  scale  1  ma  =  lv  has  again  been  used.  For  these  values 

of  L,  and  C,  ,  w  =  2 tt  ,  or  f  =  1.0  H_  .  Also,  from  the 
1  1  o  '  o  Z 

Constant-G  Chart,  for  V_  =  250  mv  and  R  =  °°,  the  amplitude 
eq  Q  o 

should  be  slightly  larger  than  250  mv  (for  sinusoidal  oscil¬ 
lation)  . 

The  curves  in  Fig.  7.4  are  the  steady-state  time- 

domain  solutions  for  V  ,  corresponding  to  the  values  of 

Gj/L./C,  as  indicated, 
dll 

It  can  be  seen  that  for  G^/L^/C^  =  *175,  the  wave¬ 
form  is  very  nearly  sinusoidal,  with  frequency  of  about  1.0 
Hz,  and  amplitude  of  about  260  mv .  These  values  are  in  good 
agreement  with  the  describing-function  theory . 

It  is  also  seen  that  for  G^/L^/C^  ~  1.0,  quite  a 
large  degree  of  distortion  is  present,  and  the  frequency 
of  oscillation  is  considerably  less  than  the  resonant  fre¬ 
quency  of  the  tank  circuit. 


v 


175  S  !  1  sec./divj  horiz 


/ 
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Figure  7 . 4 

Waveforms  Using  an  Infinite-Q  Tank  Circuit 


»  v  : 
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For  G^/l^/C-j  £  10.0,  the  solution  approaches  a 
relaxation  mode.  This  is  similar  to  the  solution  obtained 
when  the  tank  circuit  is  replaced  by  a  pure  inductor.  (See 
Ref .  3 ,  page  50 ) . 

Note  that  in  Fig.  7.4  the  peak-to-peak  amplitude 
appears  to  be  nearly  independent  of  G^/L^/C^. 


7.4  Distortion  with  Finite-Q  Circuits 


It  has  been  shown  (in  Sec.  5.9)  that  under  certain 
conditions,  increasing  by  increasing  /L^/C^  can  be  ex¬ 
pected  to  increase  distortion  and  frequency  departure.  This 
property  shall  now  be  demonstrated  by  means  of  the  two  follow¬ 
ing  examples.  Again,  in  both  cases,  the  biasing  network  has 
been  neglected,  or  =  0. 

Example  7.1 

For  this  example,  =  10,  G^/L^/C^  ~  1*0/  anc^ 
a)  =  2 tt  /l  -  1/100;  or  fQ  -  .  995  Hz.  Again  choosing 
the  current  scale  1  ma  =  1  v,  for  the  computer  simu¬ 
lation,  the  circuit  parameters  are  as  follows: 


1 

70 

1 


14  IT 


7 

2tt 


10 


7 
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Using  these  values,  the  curves  shown  in  Fig.  7.5 
(corresponding  to  the  values  of  V_.  as  indicated) 

JD 

have  been  obtained.  These  curves  represent  the 
steady-state  solutions  to  Vp.  It  is  seen  that  the 
nature  of  the  solution  depends  on  the  value  of  Vg 
chosen.  Note  that  for  all  curves  considerable 
distortion  and  frequency  departure  are  present. 

Fig.  7.5(c)  can  be  compared  to  Fig.  7.4(b)  and  (c) . 

Example  7 . 2 

In  this  case,  /l^/C^  has  been  doubled  (as  compared 

to  Example  7.1);  and  both  and  1//L^C^  have  been 

unchanged.  Thus,  =  20,  -  2.0,  and 

oo  =  2tt/1-1/400  ;  or  f  -  .999  Hz.  For  a  current 
o  '  o 

scale  of  1  ma  =  1  v,  the  circuit  parameters  for  the 
computer  simulation  will  be 

1 

70 
1 

/ 

7  TT 

7 

—  ,  and 

4  7T 

40 
7 

Note  that  R  has  been  multiplied  by  a  factor  of 
o 

four,  as  compared  to  Example  7.1. 


' 


e  J  :io^  ^xcroxxo  o  ' 


< 
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Figure  7 . 5 

Waveforms  Using  a  Tank  Circuit 
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Figure  7.6 

Waveforms  Using  a  Tank  Circuit  with  a  Q  of  20 
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Using  the  above  values,  the  solution  to  VD  (corres¬ 
ponding  to  various  values  of  V_.)  are  as  shown  in 
Fig.  7.6.  Curves  (a),  (b) ,  and  (c)  of  Fig's.  7.5 

and  7.6  can  be  directly  compared,  since  the  values 
of  VB  correspond.  The  curves  of  Fig.  7.6  clearly 
exhibit  more  distortion  and  a  higher  degree  of  fre¬ 
quency  departure  than  the  corresponding  curves  of 
Fig .  7.5. 

Curves  (d) ,  of  Fig's.  7.5  and  7.6  show  that  as  the 
peak-to-peak  amplitude  is  decreased  (by  increasing 
V  )  the  solution  more  closely  approximates  a  sinu- 

J3 

soid.  This  is  primarily  due  to  the  following  two 
factors : 

1)  as  V-.  is  increased,  r,  increases,  so  that  the 

B  d 

th 

"open-loop  gain"  to  the  n  harmonic  is  de¬ 
creased  (See  Sec.  5.9); 

2)  as  the  amplitude  of  oscillation  decreases,  the 
harmonic  content  of  the  current  wave-form  also 
decreases . 

7 . 5  Oscillator  with  Mutual  Coupling 

The  circuit  discussed  in  Example  6.1  shall  now  be 

simulated  and  tested  using  the  schematic  shown  in  Fig.  A2.16. 

For  this  circuit,  (R^  +  Rg)  =  125ft,  Q-^  =  /70,  Q ^  =  50, 

F^  =  /I . 4 ,  R  =  1000/1. 8ft,  and  V  = 
f  o  a 


237.5  mv . 


. 


' 
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Choosing 


Lici 


=  4 


TT 


(7.6) 


w  =  2tt/1  -  1/70  ;  or  f  *  .  993Hz 

°  o 


(7.7) 


From  (3.4)  and  (3.7) 


R 


o 


Q- 


=  /L1/C1  ; 


or  /L1/C1 


1000 


1.8/70 


=  66.4ft  (7.8) 


From  (7.6)  and  (7.8)  , 


66.4 


2  TT 


and 


2 TT  x  66.4 


Since 


Q-.  =  —  /C  ,  R-,  =  —  /L./C 

R-L  Q1 


66.4 


/7  0 


or  =  7.93ft.  Hence,  Rg  =  125  -  7.93  =  117.07ft. 

Again  choosing  1/ooCg  =  .01  R  ,  CB  =  lOOC^/Q^  =  11.96  C^. 

Note  that  in  Eq ' s .  A2.14,  the  only  two  quantities 
depending  on  the  secondary  tank  circuit  are  l/L^C^  and  R^/L^ 
From  (6.21)  , 


L2C2 

L1C1 


1 


or 


L2C2 


Ff  L1C1 


Thus , 


4  TT  ' 


L2C2 


1.4 


.  j\sH  bn.r»  :D?JV  I  >"  b  >  >  ^Bbroo  5  rii  no  pni  bn^c  o 
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Also , 


R, 


/L2C2  °2 


/L2C2 


r2  /c2/l2 


from  which 


R. 


2  7T 


/1.4  50 


.  04  7T 


/1.4 


For  a  current  scale  of  1  ma  =  1  v,  the  circuit 
parameter  for  the  computer  simulation  will  be: 


R±  =  .00793  , 

Lx  =  .01057  , 

C1  =  2.395  , 

Rd  =  .11707  , 

±5 

CB  =  28.6 

R  =  .556  , 

o  ' 

=  28.2  ,  and 

L2 

R 

—  =  .1062  . 


Using  the  above  values,  the  tracings  shown  in 

Fig's.  7.7  to  7.10  have  been  obtained.  Fig.  7.7  shows  the 

transient  response  of  H^(S)  (Fig.  A2.16)  to  a  unit  step  of 

voltage.  This  is  equivalent  to  applying  a  1  ma  current  step 

2 

to  the  circuit  shown  in  Fig.  6.1.  The  values  of  K  are  m- 

2 

dicated  on  the  curves.  K  has  been  adjusted  as  shown  m 
Fig.  A2.17  . 


.  saoi. 


I  <  f ■  i  'C. I  v i aqi*  :v.j  ^rislfivxups  aij  airfT 
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Figure  7 . 7 

Transient  Response  of  Mutually  Coupled  Circuits  to  a 

Current  Step 
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Voltage  Waveforms  with  Mutually  Coupled  Tank 

Circuits 


sec./div.  Horiz.  500  mv./div.  Vert.  K2  shown  on  curves 
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Figure  7 . 9 

Continuous  Solution  with  Increasing  Coupling 


: ,  :  a :  1'  S  91 


Figure  7.10 

Continuous  Solution  with  Decreasing  Coughing 
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Figure  7.8  shows  the  steady-state  solutions  to 

2  2 

VD  for  the  values  of  K  as  indicated.  Note  that  K  has  first 

been  increased  in  increments  until  a  switch  in  modes  occurred, 

and  then  similarily  decreased  until  a  second  switch  occurred. 

2 

Each  time  K  was  varied,  sufficient  time  was  allowed  for  the 

resulting  transients  to  die  out  before  the  recording  was  made. 

2 

It  was  found  that  in  varying  K  ,  the  transients 

that  were  introduced  had  a  tendency  to  precipitate  the  switch 

to  either  mode,  thus  making  it  very  difficult  to  determine 

the  exact  switching  points.  However,  the  switch  from  State 

2 

(1)  to  State  (3)  was  seen  to  occur  for  K  somewhat  larger 

2 

than  0.286,  and  the  switch  back  to  State  (1),  for  K  some¬ 
what  less  than  0.0144.  In  Example  6.1,  the  switching  points 
were  calculated  to  be  =  0.37  and  =  0.011. 

Table  7.2  compares  the  experimental  and  theoretical 

2 

values  of  amplitude  and  frequency  for  three  values  of  K  . 

The  experimental  values  have  been  taken  from  Fig.  7.8.  For 

K  =0,  the  frequency  is  given  by  (7.7);  and  the  amplitude 

has  been  determined  in  Example  6.1.  For  the  other  two  values 

of  K  ,  amplitude  and  frequency  have  been  determined  as  out- 

2 

lined  in  Example  6.1.  Note  that  K  =  0.286  corresponds  to 

2 

State  (1),  and  K  =  0.0144  corresponds  to  State  (3). 


' 


: 
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Table  7 . 2 


Comparison  of  Results 

K2 

Experimental 

A  (mv)  f(Hz) 

Theoretical 

A (mv)  f (Hz) 

0.0 

175 

.989 

177.5 

.993 

0.286 

125 

1.35 

125 

1.34 

0.0144 

125 

.834 

127.5 

.833 

In 

obtaining 

the  curves 

shown  in 

Fig ' s .  7.9 

and 

7.10,  the  pot 

2  2 

K  /1-K  , 

of  Fig.  A2 . 

17,  was 

replaced  by 

a 

servo-multiplier  with 

a  very  slow 

ramp  input  as  shown 

in 

Fig .  7.11. 


Figure  7.11 

Servo-Multiplier  Used  to  Vary  Coupling 

For  Fig.  7.9,  the  ramp  input,  x,  was  chosen  as  x(t)  =  O.Olt, 

2  2 

which  corresponds  to  a  linear  variation  of  K  /l— K  from  0.0 

2 

to  1.0,  in  104  seconds.  Thus,  for  any  one  cycle,  K  can  be 
considered  as  nearly  constant. 
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Figure  7.9  shows  the  continuous  solution  to  V^, 

9 

during  the  time  in  which  K  was  varied  from  0.275  to  0.30. 

A  few  cycles  were  missed  each  time  the  recorder  was  re¬ 
positioned  . 

Figure  7.10  is  similar  to  Fig.  7.9,  except  that  the 

system  was  initially  oscillating  in  State  (3) ,  and  K2  was 

gradually  decreased.  In  this  case,  x(t)  was  of  the  form 

x(t)  =  xo  -  O.OOlt.  Figure  7.10  shows  the  solution  to 
.  2 

during  the  time  K  was  decreased  from  0.0144  to  0.0108. 

From  Fig.  7.9,  the  switch  from  State  (1)  to  State 
(3)  appears  to  start  for  K  =  0.293;  and  from  Fig.  7.10,  the 
switch  back  to  State  (1)  appears  to  start  for  K2  -  0.0121. 

In  both  cases,  the  minimum  stable  amplitude  before  the  switch 
is  about  125  mv. 

7 . 6  Possible  Sources  of  Error 

Fig's.  7.9  and  7.10  indicate  that  the  switch  in 
modes  (from  either  state)  seems  to  occur  slightly  before  the 
theoretical  switching  points  are  reached.  This  discrepancy 
could  be  the  result  of  an  error  of  10  to  15  mv  in  VD .  Con- 
sidering  the  number  of  amplifiers  used  in  the  system,  an  er¬ 
ror  of  this  magnitude  could  be  accumulated  from  the  individual 


d.c.  offsets. 


.aJseilo  .3.6 
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It  is  also  possible  that  the  switch  from  mode  to 
mode  could  be  precipitated  by  the  presence  of  noise  in  the 
system.  As  well,  noise  at  the  input  of  the  function  genera¬ 
tor  could  effectively  alter  the  d.c.  characteristic  to  some 
extent . 


'  - 


■ 


i 
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CONCLUSIONS 

Both  singly  and  doubly  tuned  tunnel  diode  oscillators 
have  been  analyzed  in  considerable  detail,  primarily  by  means 
of  the  describing  -  function  method.  This  analysis  has  re¬ 
sulted  in  an  explanation  of  the  mode-switching  phenomenon 
exhibited  by  the  doubly  tuned  oscillator . 

On  the  basis  of  the  describing-f unction  analysis, 
a  system  was  designed,  and  then  simulated  by  means  of  an 
analogue  computer.  The  switching  phenomenon  was  indeed 
exhibited . 

The  describing-function  analysis  of  a  tunnel  diode 
oscillator  (singly  or  doubly  tuned)  is  complicated  by  the 
fact  that  the  operating  point  of  the  tunnel  diode  depends 
on  the  amplitude  of  oscillation,  and  thus  cannot  be  arbit¬ 
rarily  specified.  For  this  reason,  as  well  as  for  the 
sake  of  generality,  the  "Biasing  Curves"  and  "Constant-G 
Chart"  have  been  introduced.  These  charts  depend  only  on 
the  tunnel  diode  characteristic,  and  can  thus  be  used  in 
the  analysis  of  any  circuit  for  which  the  describing- 

function  method  is  applicable. 

A  similar  approach  could  be  taken  in  the 
describing— function  analysis  of  most  types  of  feedback 
oscillators  using  a  single-valued  non-linearity. 


. 
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APPENDIX  I 

MATHEMATICAL  JUSTIFICATIONS 


A1 . 1  Proof  of  Theorem  5.1 


The  radius  of  curvature, p,  for  a  function  y  =  f(x), 


is  normally  defined  as 

dS 


[l  +  f  (X)2J  3/2 


p 


de 


f ' ' (x) 


;  f '  '  (x)  ?  0  (Al.l) 


(Ref.  4,  page  430).  This  is  shown  in  Fig.  Al.l. 


It  will,  however,  be  more  convenient  to  consider  p  as  an 
algebraic  quantity,  whose  sign  depends  on  the  sign  of 
f ' ' (V) .  Thus,  for  the  curve  I  =  f (V) ,  p  shall  be  defined 

[l  +  f '  (V)2]3/2 


p 


f  "  (V) 


(A1.2) 
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The  quantity  dp/dS  must  now  be  found.  Since 


dp 


dS 


dp  dV 


(A1.3) 


dV  dS 


the  first  step  is  to  find  dp/dV.  Differentiating  Eq.  Al . 2 , 
dp  3f  '  '  (V)  2f  '  (V)  Ll+f  '  (V)  HI  35  ”  [l.+f  '  (V)  2]  3/2f  '  '  '  (V) 


(Al . 4 ) 


dV 

Now  since  dS^  =  dV^  +  dl^ 


1  +  f '  (V) 

dV 

dS 


f  '  '  (V) 


dV‘ 


dS  = 


1  +  f '  (V) 


or 


TT5 


dV, 


(Al . 5 ) 


[1  +  f '  (V) 

Substituting  (Al.5)  and  (A1.4)  into  (A1.3), 

dp  3f  '  '  (V)  2f  '  (V)  -  [l  +  f  '  (V)  2"]f  '  '  '  (V) 


(Al . 6 ) 


dS 


f  '  '  (V) 


This  result  has  been  cited  as  Eq.  5.9. 

Equation  Al . 6  must  now  be  evaluated  at  the  point 

( n ) 

V  =  V  .  From  Eq .  4.5,  it  is  seen  that  fy  ; (VQ)  =  n!  gR ( V^) 
Thus,  it  follows  that: 


1) 


2) 


3) 


f'(VQ)  =  gi 
f ' ' (VQ)  =  2g2 

f"'(VQ)  =  6g3 


and 


( Al . 7 ) 


Substituting  (A1.7)  into  (Al.6), 


dp 

dS 


^^1^2^  ”  3(1  +  ^1^)93 


( Al . 8 ) 


V=V 


2g 


Q 


■ 


;  ,,  •%«■>  « 


- 
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if  Vq  is  now  restricted  to  the  "negative  resistance" 
region  g^  must  be  negative.  Thus,  the  first  term  in  the  num¬ 
erator  of  (A1.8)  must  be  negative.  Furthermore,  the  denomina¬ 
tor  of  (A1.8)  must  always  be  positive  or  zero.  Therefore, 
it  is  apparent  that  if: 


1) 


dp 

dS 


^  0  , 


v=v. 


(A1.9) 


2) 


dp 


dS 


3  |  g  -^  |  j  g  3  ~ 


V=V, 


dp 

3)  — 

dS 


<  — 3 | g 1 | ,  g3  >  0 


V=V 


Q 


(A1.10) 


(Al.ll) 


If  the  range  {X}  is  now  defined  as  the  range 
Vp  <  Vq  <  Vv  ,  then  for  all  values  of  belonging  to  {X}, 
gi (Vq)  <  0  .  Assuming  that  g^ (Vq)  is  a  continuous  and  single¬ 
valued  function  in  {X},  (vq)  can  be  represented  as  shown 
in  Fig.  Al . 2 . 


range 

dp 


dS 


v=v 


Q 


The  range  {R}  can  similarily  be  defined  as  the 
<.  V  <  V-  ,  included  by  {X},  for  which 

>_  0  (if  such  a  range  exists)  .  This  is  also  illus¬ 


trated  in  Fig.  Al . 2 ,  assuming  that  dp/dS  is  a  continuous  and 
single-valued  function  in  the  range  of  interest. 

From  (A1.9),  it  is  now  apparent  that  everywhere  in 
{R},  g3(VQ)  must  be  negative  (and  non-zero). 


■ 


h':  .  sul  a  ■  .*■  "o  5  i  31  \c,5  pr.  lj-  %  £  ♦  -'  •  P *•  1  ■  - 

. 


! 


-  142 


If  it  is  now  assumed  that  g0 (V^)  is  continuous  and 

3  Q 

single-valued,  it  follows  that  there  must  be  a  larger  range, 
{ R * } ,  which  includes  {R},  for  which  g^ (Vq)  <.  0.  This  range 
is  designated  as  the  range  —  vq  —  V4 '  an<3  as  s^own  in 

Fig.  A1.2. 

It  also  follows  from  (A1.10)  that  at  , 

dp/dS  =  -3 | g1 (V^) | ;  and  at  ,  dp/dS  =  -3 | g^ (V^ ) | . 

This  proves  Theorem  5.1. 


Graphical  Representation  of 
gi(VQ),  dp/dS,  and  g3 (VQ) 


i 


.p«»*  axrfTp.O  x  ^>EP  HdWw  »S  A*  !  ■ 
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Al . 2  Elliptic  Contours 


For  a  small  region  immediately  surrounding  the  in¬ 
flection  point  of  the  curve  I  =  f (V) ,  a  reasonable  approxima¬ 
tion  to  the  function  can  be  obtained  by  using  a  cubic  equation 
of  the  form 


fT(V)  =  h  '  +  h-^V  +  h2'V2  +  h3'V3  . 
This  is  shown  in  Fig.  Al . 3 . 


Approximation  of  f (V)  by  a  Cubic  Equation 


(A1.12) 


Since  f  (V)  should  correspond  to  f (V)  at  V  =  ,  it  follows 

that : 


W 


V  (vi> 


V'V 


(VI> 


=  f(Vz)  =  go(VI)  ; 

=  f(Vz)  =  gx (Vx)  =  -Gd  ; 
=  f '  ’  (V  }  =0  ;  and 
=  f "  '  (VI)  =  6g3(VJ.)  . 


(Al .13) 
(A1.14) 
(Al .15) 


I  I  V 


(A1.16) 


. 


/  9*L/p.  i 

noidfii-pa  oxdtO  6  \d  (V>*  So  noiJemixoaqqA  1 


' 
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But  from  (A1 . 12) : 


fT'  (vz)  =  hj^'  +  2h2'VI  +  Shj'Vj  ; 


(A1.17) 


fT'  '  (Vj.)  =  2h2'  +  6h3'V].  ; 


f T  "  ’  <VI >  =  6h3'  . 


and 


(A1.18) 

(A1.19) 


Equating  (A1.19)  and  (A1.16), 


h 


3 1  =  g3(vI). 


(A1.20) 


Substituting  h3  '  into  (Al.18),  and  equating  the  result  to 
(A1.15) ; 


2h2 '  +  6VIg3(VI)  =  0  ;  or 


h2 '  = 


(Al. 21) 


Substituting  h2 '  and  '  into  (Al.17),  and  equating  to  (A1.14); 


-G 


hl'  +  2V 


h1'  -  3VJ  g3 ( VJ )  ; 


-3V];g3(VI)  +  3V].  g3(VI) 

from  which 


hi 


’Gd  +  3VI  g3(VIJ 


(Al .22) 


If  fT(V)  is  now  expanded  about  the  general  bias 
point  Vq ,  as  was  previously  done  with  f (V) ,  the  result  will 
be 


WT(V) 


fT(VQ  +  V) 


=  g 


1  +  g  'v  +  g2'v2  +  g3'v3;  (ai.23) 


. 
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where 


I  __ 


f  (n)  (V  ) 
T  K  QJ 


n 


n ! 


=  V'V 


(A1.24) 


Note  that  (A1.23)  and  (A1.24)  are  identical  to  (4.4)  and 
(4.5) ,  if  K  =  3. 

From  (4.36),  the  describing-function  in  terms  of 
the  three  coefficients  q  1 ,  will  now  be 


G  '  ( A ,  V  ) 
eq  Q' 


g,'  +  |  g3'A2 


(A1.25) 


But  from  (Al.24),  g1 '  =  f T ' (V  ) 


hl '  +  2h2'VQ  +  3h3'VQ 


(A1.26) 


Substituting  (Al.22),  (Al.21),  and  (Al.20)  into  (Al.26); 


I  _ 


-Gd  +  3VX  c,3(VT) 


+  2 


-3VIg3(VI) 


V. 


+  3g3  (V-j)  vQ 

■Gd  +  3g3(VI)  (VI2  -  2VIVQ  +  VQ2) 


•Gd  +  353<vi)(vq  -  V 


(Al .21) 


Also,  from  (Al.24), 

fm ' ' ' (V„)  6h~ ' 


I  _ 


T 


Q  _  3  _ 


=  h3'  =  g ^ ( Vj ) ;  so  that 


T  g 


'  A“ 


|  g3(vI)A2 


( Al .28) 


V 


Substituting  (A1.28)  and  (A1.27)  into  (A1.25), 


G 

eg 


(A,Vq) 


-Gd  +  3g3  (Vz )  (VQ  -  Vl)2 
+  |  g3(VI)A2  .  (A1.29) 


The  contour  G  ' (A,Vq)  =  r  is  now  defined  by  the 

equation 


•Gd  +  393(^1  (V  -  Vz)2  +  |  g3  (V-j.)  a2 


=  r ,  or 


r  +  g 


|g3(VI)A2  +  3g3  (Vj.)  (Vq-Vj)  2.  (A1.30) 


Equation  A1.30  can  be  rewritten  as 


A4 


(VQ  -  V,) 


=  1  .  (A1.31) 


(F  +  Gd> 


3g3(vI) 


(r  +  Gd) 


3g3(vi) 


Equation  A1.31  defines  an  ellipse  in  the  AVq“ 


plane,  with  centre  at  (0,VT),  major  axis’of 


r  +  G 


3g3(vi) 


and  minor  axis  of  2 


T  +  G 


3g3 (VI> 


The  major 


axis  is  parallel  to  the  A-axis.  The  above  results  have  been 
stated  in  Sec.  5.4. 

As  a  point  of  interest,  is  should  be  noted  that  if 
the  d.c.  characteristic  were  in  fact  a  cubic  (for  all  values 
of  V),  all  contours  would  be  elliptic. 


' 


- 


.9**- XXX*’  *d  M::iw  5  roinoo  IX*  '♦  ('•’  :  : 
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Al . 3  Gradient  of  G  (A,V  ) 

eg  Q' 


By 


introducing  the  notion  of  the  gradient  of  G 

eg 


v  G  ,  Eg.  5.20  can  be  derived. 

ecj 

The  differential  vector  operator,  V,  is  defined 
the  AV^-plane  as 


or 


in 


3  8 

V  =  I  -  +  j  -  ;  (A1.32) 

3  A  3V 


where  l  and  j  are  unit  vectors  parallel  to  the  A-  and  V^- 

axes  respectively.  Then  V  G  ,  or  Grad  G  ,  is  defined  as 

2  eg  eg 


V  G 


eg 


3G 

->■  eg  , 

i  - *  +  j 

3  A 


3G 


(A1.33) 


3  V 


Q 


At  some  point  P^,  V  G  can  be  graphically  inter¬ 
preted  as  a  vector,  in  the  direction  of  increasing  G  ,  and 
perpendicular  to  the  contour  G  =  r ^ ,  which  passes  through 
the  point  P^.  This  is  illustrated  in  Fig.  Al . 4 . 


-> 


Figure  Al . 4 

Graphical  Interpretation  of  the  Gradient  Vector 


' 


(  V  ■)  *o  3n:-ibeiO  C.IA 
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In  moving  the  incremental  distance,  dS  from  to 
P2'  Geq  c^an9es  from  T-,  to  T^.  This  incremental  change  will 
be  given  by 

d  G  =  V  G  *  dS  ;  (Al.34) 

eq  eq 

where 

dS  =  1  dA  +  j  dVQ  . 


Normally,  however,  the  two  partial  derivatives  of  G 
not  known,  so  that  a  graphical  means  of  determining  V 
desired . 


are 

G 

eq 


is 


If 


dS 


is  first  taken  to  be  parallel  to 


-> 

V 


G  , 

eq' 


then 


r2  -  rx  ;  or 


If  dS  is  now  replaced  by  A*S ,  which  is  a  small  but  finite 
distance  (parallel  to  V  G  ),  then  ( r ^  “  r^)  will  also  be 
small  but  finite.  Then 


|VG  I  *  — - —  (A1.35) 

e<3  I  AS  | 

On  the  Constant-G^g  Chart,  | AS |  can  be  taken  as  the  perpen¬ 
dicular  distance,  through  the  point  P^,  between  two  adjacent 

contours  G  =  rn  and  G  =  r0. 

eq  1  eq  2 


■ 
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—r 

if  VG  makes  angle  a  with  the  horizontal  axis,  it 
eg 


is  now  apparent  that 


VG 


eg 


ro  - 

2  1  ,  -t  -t-  x 

-  (l  cos  a  +  3  sm  a) 

!  AS  ! 


(A1.36) 


Since  AG  is  now  known,  the  incremental  change,  dG  ,  cor- 
eq  ^  eg 


responding  to  any  arbitrary  displacement  dS  can  be  found  from 
(A1.34) . 

If  dS  makes  angle  3  with  the  horizontal  axis,  then 


.  -t 


dS  =  i  dA  +  3  dV^  =  i  dA  +  3  tan  3  dA  =  (1+3  tan  3)dA. 


Hence,  dG  =  VG  *dS 
eg  eg 


r2  -  ri  ^ 

-  (1  cos  a  +  3  sin  a) (1  +  3  tan  3)dA,  or 


AS 


dG 


eg 


r2  -  F1 

I  AS  I 


(cos  a  +  sin  a  tan  3)dA  .  (A1.37) 


Equation  Al.37  can  be  rewritten  as 

dG 


e^ 


r9  -  r  cos(a-3) 

— - -  ( - ) 


(A1.38) 


dA 


AS 


cos  3 


In  Sec.  5.7,  dS  has  been  chosen  to  be  along  the  curve  (A) 
The  results  derived  in  this  section  can  also  be 
used  to  prove  Theorems  5.2  and  5.3. 


A1 . 4  Proof  of  Theorem  5 . 2 

Proof  of  this  theorem  is  based  on  the  gradient  of 

G  ,  as  discussed  in  Sec.  A1 . 3 . 
eg' 
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Consider  any  point,  P^  in  the  AVQ-plane.  If  dS 
is  chosen  to  be  parallel  to  the  A-axis,  then  cfs  =  1  dA,  so 
that  dV^  =  0,  and  3=0  (in  Eq.  A1.38) .  Now,  from  (A1.34) , 
dGeq  =  VGeg*dS;  and  from  (A1.33), 

^  3G  9G 

vG  =  1  — ga  +  t  gq 

GO  ^ 

9A  9Vq 

Thus,  for  dV  =0,  dG  =  9G  /9A  •  dA  or 

Q  eq  eq' 


dG  9G 

eq  eq 

dA  9  A 


(A1.39) 


axis , 
9G 

_ gq 

9  A 


If  the  point 
then  at  P^,  A  =  0. 

=  0.  Thus,  from 

A=0 


is  now  chosen  to  lie  on  the  V^- 
However,  Eq.  5.4  shows  that 

(A1.39) ,  at  P, ,  dG  /dA  =  0.  But 

1  eq' 


since  3=0,  from  (A1.38), 


dG  -  r, 

— ^  -  — - -  cos  a.  ( Al  .40) 

dA  I  AS  I 


Equation  Al.40  implies  that  cos  a  =  0,  or  a  =  ±90°.  Thus, 
all  contours  must  meet  the  axis  at  right  angles,  as 
VG  is  always  perpendicular  to  the  contour  through  the 
point  P^. 

Also,  since  moving  in  any  direction  from  the  point 
( 0 , )  corresponds  to  increasing  G  ,  it  follows  that  if 
V  >  Vj ,  a  =  +90°,  and  if  <  V^. ,  a  =  -90°.  This  is  shown 
in  Fig  .  Al  .5  . 


w-  - 


j  i.  cq  riJ  /  t  i  n  ..b  Y  6  Pr  ,/rj  1  °  ' 
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A1 . 5  Proof  of  Theorem  5.3 


Proof  of  this  theorem  is  based  on  the  results  ob¬ 
tained  in  both  Sec's.  Al . 3  and  A1 . 4 . 


82G 


eq 


3  A 


Consider  a  curve  G  (A,V_  ) ,  such  that 

eg  Q0 


<  0.  This  curve  will  clearly  concave  downward  at 


A=0 


A  =  0.  One  such  curve  is  shown  in  Fig.  5.4(a)  ,  and  again 

3G 

eq 


in  Fig.  Al . 6 .  Now  since 


3  A 


=  0  (Eq.  5.4) ,  it  follows 


A=0 


that  there  must  exist  a  small  positive  number  e,  for  which 


3G 


eq 


3  A 


<  0 


A=e 


'  '  .  .  ;o« 
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If  the  point  (e,Vn  )  is  now  chosen  as  in  the  AV 


and  dS  is  chosen  parallel  to  the  A-axis,  then,  in  Eq 

dG  3G 

6  =  0,  and  by  (A1.39)  dAe  -  = 

A=e 


plane , 
Al. 38, 
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Thus,  at  ,  dG  /dA  <  0.  But  since  3=0, 


dG 


dA 


r2  -  h 

i  AS  I 


cos  a,  which  implies 


that  a  >  +90°,  or  a  >  -90°. 

If  it  is  now  assumed  that  V  >  V  ,  it  is  apparent 

g2  1 

that  by  making  e  arbitrarily  close  to  zero,  a  can  be  made 
arbitrarily  close  to  +90°.  (See  Fig.  A1.5).  Hence,  it  can 
be  concluded  that  at  P^,  a  >  +90°,  which  implies  that  the 
contour  through  must  have  a  positive  slope  at  P^.  This 
is  shown  in  Fig.  A1.7. 

Since  it  is  known  that  all  contours  meet  the  V^-axis 
at  right  angles,  and  since  e  can  be  made  arbitrarily  close 
to  zero,  it  now  follows  that  the  contour  through  P^  must  be 
concave  upward  as  A  approaches  zero. 

However,  it  can  also  be  seen  that  this  contour  will 
not  meet  the  V  -axis  at  V  ,  but  at  a  slightly  smaller  value. 

y  y2 

If  this  value  of  V  is  defined  as  (V  -  £) ,  it  is  apparent 

y  y2 

that  as  e  is  made  arbitrarily  small  (but  finite)  the  same 
will  be  true  of  £.  For  a  specific  e,  £  will  of  course  de¬ 
pend  on  the  degree  of  curvature  of  the  contour  in  the  vicin¬ 
ity  of  P-^. 

If  the  range  {V  }  is  now  defined  as  the  range 

U2 

V.,  <_  V  <_  V9 ,  such  that  Vi  >  V  ,  and  such  that  all  curves 

1  Q  2  II 

G  (A,V^)  corresponding  to  {V^  }  are  concave  downward  at 
eq  Q  ^  Q2 

A  =  0,  then  it  can  be  concluded  that  all  contours  meeting 


. 


9C,&e  (si -an  i,<  t  U*£*  fits  sb6m  8  :  “  : 

' 


' 
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the  V^-axis  in  the  range  -  ^2.  ~  VQ  “  V2  “  ^2  must  ke  con_ 

cave  upward  at  A  =  0.  Now  by  making  e  arbitrarily  small, 

and  £ 2  will  also  be  made  aribtrarily  small,  so  that  the 

range  <_  <_  V2  -  £2  can  be  made  to  correspond,  as 

closely  as  desired,  to  the  range  {V^  }.  The  end  points  of 

Q2 

the  two  ranges,  however,  cannot  be  made  to  coincide,  since 
e,  and  thus  and  £2  must  remain  finite. 

A  similar  argument  can  be  used  to  prove  the  re¬ 
mainder  of  Theorem  5.3. 


-nos  ed  ieum  .}  -  V  t  7  z  i  -  ,V  spn«:s  arid  .x  ;  txf.-0'  aria 

arid  3 o?  ,  X I  ama  v  E i  *  4*3 *d  x  i /:  obis  tt  sd 


’ 
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APPENDIX  II 

ANALOGUE  COMPUTER  SIMULATION 

It  has  been  shown  in  Sec.  7.1,  that  the  tunnel  diode 
oscillator  (with  or  without  mutual  coupling)  can  be  accurately 
simulated  on  an  analogue  computer,  by  means  of  three  basic 
system  components.  Computer  schematics  for  these  basic  com¬ 
ponents  will  now  be  derived. 

A2 . 1  Function  Generator 

A  function  generator  which  will  produce  the  d.c. 
characteristic,  I  =  f (V) ,  will  first  be  designed.  This  is 
done  by  using  an  actual  tunnel  diode,  so  that  the  resulting 
function  will  be  as  close  as  possible  to  the  actual  curve 
for  the  diode.  This  function  generator  has  been  used  to  ob¬ 
tain  the  characteristic  shown  in  Fig.  4.1(a). 
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with  open 
back  resi 
the  input 

less  than 
i2  -  Ae/R 


so  that 


Figure  A2 . 1  shows  an  ordinary  operational  amplifier, 
-loop  gain  of  A,  input  impedance  of  Z^,  and  a  feed- 
;tor,  R  .  A  tunnel  diode  is  placed  in  parallel  with 
resistor ,  R . . 

i 

If  it  is  assumed  that  the  grid  voltage,  e,  is  much 
or  -V2 ,  t*1611  ~  V^/R^  +  f  (V^)  ,  and 

Since  e  =  (i^  -  i^)Z it  follows  that 


V. 


Ae 


—  +  f  (V  )  - - 

R.  R, 

l  f 


Z .  ;  or 
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Z  . 

V, 

e 

1  +  A(— ) 

a 

—  +  f  (Vn) 

Rf_ 

R.  1 

1 

Z ^  .  But  e  =  -V^/A, 


V. 


Z  . 

V1 

1  +  A(— ) 

a 

—  +  f  (V,) 

r£ 

R.  1 
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V, 


AZ  . 
1 


1  +  A(— ) 
R- 


vx 

—  +  f (V.) 

R. 

1 


If  it  is  further  assumed  that  A(Zi/Rf)  >>  1,  then 

V, 


V. 


-  R 


R 


+  F(V1) 


=  -  (—  )V1  -  Rff  (Vx)  . 


R. 

1 


If  the  quantity  (Rf/Ri)V1  is  now  added  to  V2,  the 
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result  will  be  -R^f(V^).  This  can  be  done  as  shown  in  Fig. 
A2 . 2. 


Figure  A2 . 2 

The  Function  Generator 


The  reason  for  including  ,  in  parallel  with  the 
tunnel  diode,  is  to  prevent  self-oscillation  of  the  diode 
(in  conjunction  with  the  circuitry  of  the  analogue  computer) , 
as  discussed  in  Sec.  5.10.  To  assure  that  this  does  not  oc¬ 
cur,  R^  can  be  connected  directly  to  the  leads  of  the  tunnel 
diode  (so  that  parasitic  elements  will  be  minimized);  and 
can  be  chosen  in  magnitude  so  that  the  resulting  combination 
does  not  exhibit  a  ’’negative  resistance”  region.  Thus,  from 
Fig .  A2 . 1 , 

d.  -| 

^"1  ^ 

— -  >0,  or  —  >  -f ' (V  ) . 

dV.  R. 

l  i 


. 
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Since  the  maximum  value  of  -f ' (V^)  is  (Eq.  1.2), 

this  implies  that  1/R.  >  G, ,  or  R.  <  1/G  ^  . 

id  1  '  d 

For  a  1N2940  tunnel  diode,  G^  is  typically 

-3 

6.6  x  10  mho  (Table  4.1)  ,  so  that  R^  should  be  less  than 
150  n.  R^  can  be  chosen  as  any  desirable  value. 


A2 . 2  Simulation  of  Biasing  Network 


The  transfer  function,  (S)  ,  mathematically  equal 

to  the  biasing  network,  Z_  ( S )  ,  will  now  be  simulated.  (This 

a 

is  discussed  in  Sec.  7.1) 

From  Fig .  2.1  (b)  , 


Vs)  = 


RB  ’  1//SCB 

RB  +  ^SCB 


R 


B 


scbrb  +  1 


=  h2(s) 


The  computer  schematic  for  this  transfer  function  can  very 
easily  be  derived,  and  is  shown  in  Fig.  A2 . 3 . 


Figure  A2 . 3 

Simulation  of  the  Biasing  Network 
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A2 . 3 


H1(S) 

tank 


Equat 


where 

H1(S) 


Simulation  of  Single  Tank  Circuit 

The  computer  schematic  for  the  transfer  function 
,  equal  to  the  impedance  function,  ZT(S)  for  a  single 
circuit  will  now  be  derived. 

From  Eq.  3.1, 


ZT(S)  = 


R^  +  SL^ 


S  +  SC^R-^  +  1 


=  H1(S) 


(A2.1) 


Lon  A2 . 1  can  be  rewritten  as 

R. 


Hn  (S)  =  - 

^  C 

U1 


s  + 


2  ^1  ^ 
Sz  +  _  S  +  - 


L1C1 


1 

c. 


S  +  a 


4-  b ,  S  +  b 
1  o 


(A2.2) 


a  ,  b-,  ,  and  b  are  implicitly  defined.  Letting 
o  1  o 

=  C^(S)/R^(S),  (A2.2)  can  be  rewritten  as 


c1(s) 

R-L  (S) 


1 

C. 


S  1  +  a  S  2 
o 


1  +  bxS  1  +  b1S  2 


or 


cps)  = 


I^tS) 


S  1  +  a  S  2 
_ o _ 

1  +  bn  S_1  +  b  S-2 

1  o 


(A2.3) 


If  E  ( S )  is  defined  as 


E  ( S )  = 


R±  (S) 


Cld  +  b^-1  +  boS’2) 


(A2.4) 
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(A2.3)  can  be  expressed  as 


C1(S)  =  E(S)  (S  1  +  aQS  2) 

Now,  (A2.4)  can  be  rewritten  as 


(A2.5) 


-1  9  Ri  <S) 

E(S)  (1  +  b.S  1  +  b  S  Z)  =  — - 

1  °  C 

C1 


R1(S)  -1  -? 

E  ( S )  = - b.  S  E(S)  -  b  S  E(S)  . 

r  1  o 


or 


(A2.6) 


Using  the  method  of  Signal  Flow,  (A2.6)  can  be  re 
presented  as  shown  in  Fig.  A2 . 4 . 


Similarily,  (A2.5)  can  be  represented  as  shown  in  Fig.  A2 . 5 


E(S) 


cx  (S) 


Figure  A2 . 5 

Signal  Flow  Graph  for  Equation  A2 . 5 


. 


-  161 


Combining  Fig's,  A2 . 4  and  A2 . 5  results  in  the  complete  Signal 
Flow  Graph  for  (S)  ,  which  is  shown  in  Fig.  A2 , 6 . 


-b 


Figure  A2 . 6 

Signal  Flow  Graph  for  a  Single  Tank  Circuit 

The  computer  schematic  can  now  be  drawn  directly 


from  Fig.  A2 . 6 ,  and  is  shown  in  Fig.  A2 . 7 .  The  values  of 
aQ r  /  and  b^  (as  defined  by  (A2.2))  are  shown  on  the 
schematic . 


Figure  A2 . 7 

Computer  Schematic  for  a  Single  Tank  Circuit 


Figures  A2 . 2 ,  A2 . 3  ,  and  A2 . 7  can  now  be  combined 
to  form  the  complete  simulation  of  the  circuit  shown  in 
Fig.  2.2.  The  resulting  schematic  is  shown  in  Fig.  A2 . 8 


(See  Fig .  7.1). 


.  .  '  1  ■ 


. 

' 
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Figure  A2 . 8 

Computer  Schematic  for  an  Oscillator 
Using  a  Single  Tank  Circuit 


, 
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It  was  suggested  in  Sec.  5.9,  that  the  dependence 
of  distortion  on  the  factor  G^/L^/C^  could  be  investigated  by 
simulating  a  tunnel  diode  oscillator  with  an  infinite-Q  tank 
circuit.  From  Fig.  A2 . 8 ,  this  can  easily  be  done  by  discon¬ 
necting  the  two  pots,  R^/L.,  and  R^/L^C^. 

A2 . 4  Simulation  of  Mutually  Coupled  Tank  Circuits 

If  the  impedance  function,  Z^CS),  as  given  by  Eq. 

6.2,  is  directly  synthesized,  as  was  done  for  Eq.  A2 . 1 ,  it 

2 

has  been  found  that  the  quantity  (1  -  K  )  will  appear  m 

several  different  parameters  on  the  resulting  computer 

schematic.  Hence,  a  different  approach  will  initially  be 

2 

taken,  and  it  will  be  seen  that  K  can  be  varied  by  adjust¬ 
ing  a  single  pot. 


(a)  (b) 

Figure  A2 . 9 

Circuits  with  Mutual  Coupling 

It  is  desired  to  simulate  a  transfer  function  H^S), 
equal  to  ZT(S)  =  V(S)/I(S),  as  shown  in  Fig.  A2.9(a). 


' 
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Consider  first  the  circuit  shown  in  Fig.  A2.9(b),  where 
Z2(S)  =  V2(S)/I  (S).  This  circuit  is  identical  to  that 
shown  in  Fig.  6.2,  except  that  has  been  added.  Thus, 
Z2  (S)  =  R^  +  Z-^(S),  where  Z^(S)  is  defined  by  Eq.  6.8;  or 


SL. 


Z  2  ( S )  =  R±  + 


S2L2C2 (1-K2)  +  sr2C2  + 


(S  L2C2  +  SR2C2  +  1) 


This  can  be  rewritten  as 


z2(s)  = 


S3L1L2C2  (1-K2)  +  s2  (L1 R2C7+R1L2C2)  +  S(L1+R1R2C2)  +  R-|_ 


(s  l2c2  +  sr2c2  +  1) 


(A2.7) 


Now  let 


Y2(S)  = 


Z2(S) 


N(S) 
D  ( S ) 


(A2.8) 


where  N(S)  and  D(S)  are  the  denominator  and  numerator  of 
(A2.7)  respectively.  Also,  letting  Y,^(S)  =  1/Z^,(S),  from 
Fig.  A2 . 9 (a)  it  follows  that  YT(S)  =  Y2(S)  +  SC1,  or 


yt(s)  =  sc1  + 


N(S) 
D  ( S ) 


SC^tS)  +  N(S) 
D  (S ) 


Hence , 


ZT(S) 


be  rearranged  as 


D(S) 


yt(s) 


SC1D(S)  +  N(S) 


which  can 


zT(s) 


SC. 


N(S) 


1  + 


SC1D ( S ) 
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Thus,  the  problem  is  now  to  simulate  the  transfer  function 


H1(S) 


1 

SC 


1 

N(S) 


1  + 


C1(S) 

Rx(s) 


SC1D(S) 

The  closed  loop  system  shown  in  Fig.  A2.10  will  have 
transfer  function. 


this 


Mutually  Coupled  Tank  Circuits 


The  feedback  portion  of  Fig.  A2.10  can  now  be 
simulated  by  the  same  means  as  was  Ecr.  A2 . 1 .  From  Eq.  A2 . 7 , 


N(S) 


S  L2C2  +  SR2C2  +  1 


d (s)  s3l1l2c2 (1-k2) +s2 (l1r2c2+r1l2c2) +s (l1+r1r2c2) +R 


(A2.9) 


which  can  be  written  as 


N  ( S ) 


D(S) 


2  R2  ^ 
S  +  S  =-=-  + 


L2C2 


s3(1-k2)+s2(j-^  +  Fi)+s(r-5- 

L2  L1  L2c2 


R1R2 

+  YTir)  + 

L1L2 


R 


1 


L1L2C2 


(■ 


S  +  a,  S  +  a 
1  o 


L±  b3S3  +  b2S2  +  bxS  +  bQ 


)  . 


(A2.10) 
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Letting  N(S)/D(S)  =  C(S)/R(S),  (A2.2)  can  be  expressed  as 

C(S)  1 


R(S) 


-1  -2  -3 

S  +  a,  S  +  a  S 
_  1 _ o _ 

L,  b~  +  b 0 S  1  +  b, S~2  +  b  S-3 
13  2  1  o 


(• 


)  ;  or 


R (S)  S  1  +  a1  S  2  +  a  S  3 

C  ( S )  =  -  ( -  1  - -~)  .  (A2.ll) 

L,  b0  +  b0S  +b,S+bS 
13  2  1  o 


Defining  E(S)  as 


E(S)  = 


R  (S) 


L, (b.  +  bnS  1  +  b,S  2  +  b  S  3) 
13  2  1  o 


(A2.12) 


(A2.ll)  can  be  written  as 


C(S)  =  E  ( S )  (S  1  +  a1S  2  +  aQS  3) 


(A2.13) 


Also ,  (A2.12)  can  be  rewritten  as 


E(S)  (b,  +  b2S  1  +  b1S  2  +  bQS  3)  = 


R  (S) 


or 


R  (S) 


b3E(S)  = 


-  b2S  1E(S)  -  b1S  2E(S)  -  bQS  3E(S)  (A2.14) 


Equation  A2.14  can  be  represented  by  the  Signal  Flow 


Graph  shown  in  Fig.  A2.ll 


*  ^ 
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Figure  A2.ll 

Signal  Flow  Graph  for  Ecmation  A2.14 
Also,  Eq.  2.13  can  be  represented  by  Fig.  A2.12. 


Figure  A2.12 

Signal  Flow  Graph  for  Equation  A2.13 

Combining  Fig's.  A2.ll  and  A2.12,  the  Signal 
Flow  Graph  for  the  feedback  function  N(S)/D(S)  results. 

2 

This  is  shown  in  Fig.  A2.13.  Note  that  b^,  which  is  (1-K  ), 
appears  only  once  in  this  figure. 
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Figure  A2.13 

Signal  Flow  Graph  for  Feedback  Portion  of  Fig.  A2.10 

The  computer  scheme  for  N(S)/D(S)  can  now  be 
drawn  from  Fig.  A2.13,  and  is  shown  in  Fig.  A2.14. 


Figure  A2.14 

Computer  Schematic  for  Feedback  Portion  of  Fig.  A2.10 

From  Fig.  7.1,  the  block  diagram  for  the  complete 
simulation,  with  mutual  coupling,  will  now  be  as  shown  in 


. 
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Fig.  A2.15.  Thus,  the  complete  computer  schematic  will  be  as 
shown  in  Fig.  A2.16.  (From  Fig's.  A2 , 2 ,  A2 . 3 ,  A2.14,  and 
A2.15) . 


B 


D 


Figure  A2.15 

Block  Schematic  for  Oscillator  with 
Mutually  Coupled  Tank  Circuits 


For  convenience,  the  coefficients  a  and  b  are 

n  n 


listed  below. 


L2C2 


R. 


R. 


L1L2C2 


R1R2 


L2C2  L1L2 


^2  +  ^1 

L2  L1 


(A2.14) 


b„  =  (1  -  K  ) 


■  / .  -  : 


!  nsr  ■  ■  -  ' 


m 


m  -  x) 
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Figure  A2.16 

Computer  Schematic  for  Oscillator  with 
Mutually  Coupled  Tank  Circuits 
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A2 . 5  Accurate  Adjustment  of  K‘ 


Note  that  in  Fig's.  A2.14  and  A2.16,  the  coeffici- 

2 

ent  l/b^  represents  the  quantity  1/(1-K  ),  so  that  varying 

2 

the  pot  1/10 (1-K  )  from  0.1  to  1.0  corresponds  to  varying 
2 

K  from  0.0  to  0.9.  However,  in  practice,  difficulty  would 

2 

be  encountered  in  accurately  adjusting  K"  to  very  small  values, 

2 

or  to  zero.  Thus,  an  alternate  means  of  varying  K  would  be 
desirable . 

2 

The  quantity  1/(1-K  )  can  be  expressed  as 

.2 


K 


=  1  + 


1-K  1-K 

Thus,  Fig.  A2.16  can  be  modified  as  shown  in  Fig.  A2.17, 
where  the  amplifier  and  integrator  shown  in  Fig.  A2.17  re¬ 
present  the  amplifier  and  integrator  connected  to  the  pot 


1/10 (1-K  ),  in  Fig.  A2.16 


2  2 
K  / (1-K  ) 


Figure  A2.17 

2 

Accurate  Adjustment  of  K 

Hence, varying  the  pot  K2/(l-K2)  from  0.0  to  1.0  corresponds 

9  2 
to  varying  K  from  0.0  to  0.5.  For  very  small  values  of  K  , 

K2/(l-K2)  ^  K2,  and  for  K2  =  0,  the  pot  can  be  disconnected. 


. 


. 
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A2 . 6  Scaling 

For  the  computer  simulation  of  any  practical  oscil¬ 
lator  circuit  (with  or  without  mutual  coupling)  frequency 
scaling  will  be  required.  This  can  be  done  by  simply  scal¬ 
ing  the  impedance  functions,  ZT(S)  and  Z  (S)  ,  to  any  desired 
frequency,  and  using  the  resulting  circuit  elements  to  deter¬ 
mine  the  pot-settings  and  gains  for  the  computer  schematics. 

It  is  shown  by  Skilling  (Ref.  2,  page  501)  that 
L  C 

Z(R,  -  ,  -  ,  nS)  =  Z (R , L , C , S )  (A2.15) 

n  n 

Thus,  if  it  were  desired  to  simulate  a  circuit  that  was 
known  to  oscillate  at  1.0  MHz, and  the  simulation  was  desired 
to  oscillate  at  1.0  Hz,  in  Eq.  A2.15,  n  =  10  Hence,  for 

the  simulation,  all  resistors  would  be  unchanged,  and  all  in- 

g 

ductors  and  capacitors  would  be  multiplied  by  10  . 

Impedance  scaling  can  also  be  applied  if  desired. 

If  all  linear  impedances  are  multiplied  by  a  factor,  m,  the 
form  of  the  solution  will  remain  unaltered,  provided  that  the 
diode  characteristic,  f (V) ,  is  multiplied  by  1/m. 

It  should  be  kept  in  mind  that,  on  the  computer 
simulation,  is  the  actual  voltage  appearing  across  the 
tunnel  diode,  and  thus  must  not  be  allowed  to  exceed  the 
operating  range  of  the  diode.  Any  desirable  scale  can  be 
chosen  for  X2,  which  represents  the  diode  current. 


. 
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